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^ ■ ABSTRACT 

/^s | We investigate the interplay of cosmic ray (CR) propagation and advection in galaxy clusters. Propagation in form of CR diffusion 

and streaming tends to drive the CR radial profiles towards being flat, with equal CR number density everywhere. Advection of CR 
\uJ ■ by the turbulent gas motions tends to produce centrally enhanced profiles. We assume that the CR streaming velocity is of the order 

of the sound velocity. This is motivated by plasma physical arguments. The CR streaming is then usually larger than typical advection 
velocities and becomes comparable or lower than this only for periods with trans- and super-sonic cluster turbulence. As a consequence 
a bimodality of the CR spatial distribution results. Strongly turbulent, merging clusters should have a more centrally concentrated CR 
energy density profile with respect to relaxed ones with very subsonic turbulence. This translates into a bimodality of the expected 
diffuse radio and gamma-ray emission of clusters, since more centrally concentrated CR will find higher target densities for hadronic 
CR proton interactions, higher plasma wave energy densities for CR electron and proton re-acceleration, and stronger magnetic fields. 
Thus, the observed bimodality of cluster radio halos appears to be a natural consequence of the interplay of CR transport processes, 
independent of the model of radio halo formation, be it hadronic interactions of CR protons or re-acceleration of low-energy CR 
electrons. Energy dependence of the CR propagation should lead to spectral steepening of dying radio halos. Furthermore, we show 
that the interplay of CR diffusion with advection implies first order CR re-acceleration in the pressure-stratified atmospheres of galaxy 
clusters. Finally, we argue that CR streaming could be important in turbulent cool cores of galaxy clusters since it heats preferentially 
the central gas with highest cooling rate. 
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<^> 1. Introduction I1997I) . CRe suffer much more severe energy losses via syn- 

t-H . chrotron and inverse Compton emission at GeV energies, and 

. 1. 1. Motivation Bremsstrahlung and Coulomb losses below 100 MeV. CRe with 
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Relativistic particle populations, cosmic rays (CR), are expected an ener gy of ~ 10 GeV emit GHz synchrotron waves in ^G- 
to permeate the intra-cluster medium (ICM). Cosmic ray elec- stren g th magnetic fields. Since the associated inverse Compton 



b ■ irons (CRe) are directly visible in many galaxy clusters via and synchrotron cooling time is r IC , sy „ ~ 2 x 10 8 yr, these CRe 

their radio synchrotron emission, forming th e so-call ed clus- must have been recentl y ln J ected or re-accelerated, whereas any 

ter radio halos (e.g. Kemrjne7irid]£004; FerettIi[d]SQ04b). CR P can be as old as lts cluster Slnce CR P should be ab undant 

Several CRe injection sites can also be identified via the same and are able to ln J ect relativistic electrons via the production 

synchrotron radiation mechanism: winds and gas stripping from of char 8 ed P lons ln hadronic interactions with gas nuclei, they 

cluster galaxies, active galactic nuclei, and shock waves from could be the on g ln of the observed radio halo^ Alternatively, 

structure formation. All these should also be injection sites for a low ener gy ( 10 ° MeV > CRe ICM Population might be suffi- 

CR protons (CRp) and heavier relativistic nuclei. Due to their Clentl y lon § llved lf 11 18 maintained by re-acceleration by plasma 

higher masses with respect to the electrons, protons and nuclei waves a 8 ainst coolln g Processes. During the phases of high ICM 

are accelerated more efficiently. In our own Galaxy, the ratio of turbulence after cluster merger the re-acceleration might be so 
the spectral energy flux of CRp to CRe between 1. . . 10 GeV is 
about one hundred (ISchhckeiserl I2002T) . Similar ratios are also 
expected at least for the injection from galaxies and structure 
formation shock waves for the same kinematic reasons. 

Cluster CRp should have accumulated over cosmic i Dennispn fi^g^^ fi^. | B lasi & Colafrancescol 



„ _ _ „ (1999); 

timescales since the bulk of them is unable to leave through the Do i ag & EnBlin 12000): iMiniati et all fcOOlbl): iMiniatil 1 2003) 



persistent infall of matter onto the cluster and due to the long pf romm er & EnBlin i2004al lbl): IPfrommer et al.1 d2008h: iPfrom mer 
CRp's radiative lifetime s in the ICM of the o r der of an Hubble |2008h ; |Kushnir et aljBo09l) : lDonnert et alj j2010al lbh; lKeshet & Loeb 
time (IVolk et all Il996t lEnfilin et al.1 119971: iBerezinskv et all d20ld> : lKeshetl J2010h 
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efficient that CRe are accelerated into the radio observable en- 
ergy range of ~ 10 GeVQ 

For both particle populations, CRp3 and CRe, which we ad- 
dress commonly as CR, the question is what spatial and spectral 
shape they have acquired? This will largely determine which ra- 
diative signatures and which dynamical influences we can expect 
from their presence. 

Without continuous injection of CRe, their spectra should be 
highly curved, due to the severe energy losses relativistic elec- 
trons suffer from, both at the high and low end of their spectra. If 
not replenished, or sufficiently re-accelerated, the CRe popula- 
tion should disappear in the cluster center after a time of about a 
Gyr. Thus, the CRe visible in the radio halos should be spatially 
close to their injection or last re-acceleration site. 

At energies above a few GeV, the injected CRp power law 
spectra are unaffected by hadronic and adiabat ic losses, only 
the normalisation of the spectra evolves (e. g., iMiniatil [2001; 
lEnBlin et alJ 12001 iPinzke & Pfrommeril2010l) . The spatial dis- 
tribution, however, may be strongly affected by transport pro- 
cesses. Also local spectral modifications can be expected if the 
macroscopic transport mechanism depends on the CR energy. 
Although the CRp travel close to the speed of light, the ICM 
magnetic fields of O(fiG) force them on helical orbits along the 
field lines, and resonant scattering events with plasma waves will 
try to isotropise the CR momentum distribution rapidly. Thus 
CRp are tied to the fields and their transport is strongly con- 
trolled by the magnetic and turbulent properties of the ICM. 

In addition to the radiative signatures of CRs in the radio and 
gamma ray bands, CRs can store significant amount of energy, 
which they preferentially release in the centers of clusters due to 
the higher target density there and because of the adiabatic losses 
they suffer when they propagate to the cluster outskirts. Thus, 
CRs were proposed to help to stabilise cluster cool cores against 
a cooling instability^. The spatial distribution of the CR popu- 
lation, as shaped by transport processes, as well as the energy 
absorbed and released during advective and streaming transport 
are of direct importance for this. Three transport processes are 
relevant here: 

- Advection: The magnetic field lines are largely frozen into 
the thermal plasma of the ICM and are dragged with any 
gas flow. The enclosed CR are advected with this flow and 
suffer energy losses or gains from any adiabatic expansion or 
compression of the flow, respectively. 

- Diffusion: A CR may travel several gyro-radii along a field 
line before it is resonantly scattered by plasma waves of 
the medium. The resulting random walk along the field line 
leads to a considerable diffusion parallel to the field lines. 
Since the gyro center of the particles are displaced in the 
plasma-wave interaction, a small perpendicular diffusion re- 
sults. This leads to a larger macroscopic displacement from 
the original magnetic field line when the particles follow the 

2 ISchlickeiser et al.1 dl987h: iGiovannini et all Jl993h : ferunetti et ail 
d2004h: IBrunetti & Blasil d2005h : iBrunetti & Lazarianl d2007l 1201 Oh : 
IBrunetti et al.ld2009h 

3 Actually, there should also be a substantial population of relativis- 
tic alpha-particle. These can, however, regarded to be equivalent to an 
ensemble of four CRp, traveling together due to the relatively weak 
nuclear bin ding forces betwee n them (and keeping two of them as neu- 
trons). See EnBlin et al. (2 0071) for an extended discussion of this. 

4 Heating via Coulomb lo sses of CRs fr om a cen- 

tral AGN are discussed i n | Colafrancesco et all {2004) and 
IColafrancesco & Marchegianil d2008t) . T he additional he ating via 
streaming of the same CRs is considered in lGuo & Ohl d2008h . 



diverging path of this initially neighboring line. The CR dif- 
fusion coefficient generally increases with the particle en- 
ergy. 

- Streaming: In the presence of a sufficiently large gradient in 
the spatial CR distribution along a field line, an anisotropic 
momentum distribution function builds up, since more par- 
ticle arrive at such a location from one side than from the 
other. This leads to a net CR flux towards the CR depleted 
region. The streaming of the CR with respect to the thermal 
ICM excites plasma waves, on which the particles scatter. 
This limits the streaming velocity of CRs to be of the order 
of the Alfvenic or the sound speed, depending on the plasma 
properties. 

Advection is a passive form of transport, which is actually in- 
cluded in many of the numerical simulation schemes for CRp 
in the large-scale-structur^- Diffusion and streaming are ac- 
tive propagation processes, relying on the own speed of the CR. 
Inclusion of CR diffusion exists for some MHD codefl which 
are, however, not suited for large-scale-structure simulations. To 
our knowledge, CR streaming is not implemented in any of the 
simulation codes used in astrophysics. Thus, the effects we dis- 
cuss in this paper concern physics, which is not captured in cur- 
rent numerical simulations of galaxy clusters. 

1.2. Goal and outline 

With this paper, we want to discuss the possible effects CR 
propagation has on the spatial and spectral distribution of CR 
in galaxy clusters and their radiative signatures in gamma-rays 
and radio frequencies. At this stage of the research, without 
the necessary cosmological simulation tools including simulta- 
neously CR propagation and magnetic field evolution at hand, 
it is not possible to make definite predictions, due to the com- 
plexity of the interaction of turbulent ICM gas with CR and 
magnetic fields. Our goal here is to outline plausible scenarios, 
which highlight the potential importance of CR propagation, and 
thereby motivates further research and hopefully lead to the de- 
velopment of the necessary simulation tools. 

A strong motivation for such developments should be the 
fact that CR propagation can explain the observed bimodality of 
cluster radio halos within the hadronic radio halo model. Since 
all necessary ingredients of this model, as CRp shock accelera- 
tion, a long CRp cooling times, the processes of hadronic e ± pro- 
duction and e ± synchrotron losses, as well as the required field 
strength of a few fiG, are known to be given in typical ICM envi- 
ronments, this model is very natural, nearly free of assumptions 
and therefore attractive. The main counter-argument, that not ev- 
ery cluster exhibits a radio halo, is alleviated if CR propagation is 
operative in the ICM as argued here. However, also radio halos 
in the CRe re-acceleration model should strongly benefit from 
propagation effects. The low energy CRe population, which is 
re-accelerated to explain radio halos in these models, is prob- 
ably thermalised in the cluster centers due to Coulomb energy 
losses during quiet phases of the cluster, when re-acceleration 
is weak. However, CRe can survive for Gyrs in the outskirts of 
clusters during such phases, and being dragged into the core dur- 
ing later turbulent merger phases. Then, the required turbulent 
re-acceleration is also present, in order to power a radio halo in 



5 IMiniatil d2001l I2002L I2007T): iMiniati et all d2007l) ; iPfrommer et all 
d2006l);lEnfilinetal.ld2007h : lJubelgasetal.ld2008l) 

6 iBrandenbure & Dobleil d2002l); iBrandenburd 420031) ; 
lHanasz & Leschl d2003h : ISnodin et al] d2006h : lHanasz et al] feOlojbh 
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the re-acceleration model. Also for this model, the observed ra- 
dio bimodality of clusters should partly result from the interplay 
of the CRe transport mechanisms. 

The outline of this work is the following. In Sect. [2j we first 
introduce radio halos, their connection to the dynamical state of 
clusters, and the strength and weaknesses of the two main theo- 
retical scenarios used to explain them. Then we discuss in Sect. [3] 
typical conditions of the turbulent ICM and their implications for 
CR transport processes and the implied adiabatic energy gains 
and losses. The consequences for non-thermal cluster emission 
in the gamma-ray and radio spectral bands as well as the poten- 
tial heating of cool cores are investigated in Sect. |4] Finally we 
conclude in Sect. [5] 

2. Radio halos 

2.1. Observational properties 

Cluster radio halos are our primary evidence for the existence 
of CR in galaxy clusters. They are spatially extended regions of 
diffuse radio emission, which have regular morphologies, very 
much like the morphology of the X-ray emitting thermal ICM 
plasma. Their radio synchrotron emission is unpolarised, due to 
the contribution of various magnetic field orientations along the 
line of sight, and Faraday rotation de-polarisation. A compilation 
of radio halo luminosities and X-ray properties of their hosting 
clusters can be found in Tab. Q] and Fig.Q] 

Cluster radio halos have t o be discriminated fr om cluster ra- 
dio relics (for taxonomy see iKempner et al]|2004f) . Relics are 
also extended, but often located at the cluster periphery of galaxy 
clusters. They show an irregular morphology and are often sig- 
nificantly polarised. As halos, they appear preferentially in clus- 
ters exhibiting signs of merging activity. A halo and relic can 
appear in the same clustefl or even two relics simultaneously 
have been observecfj Relics are believed to trace merger and 
accretion shock waves, either due to direct particle acceleration 
via the Fermi-I mechanism^, or due to the revival of old radio 
cocoons of radio galaxies via compressiorFl. Multi relic for- 
mation emerges naturally in numerical models of structure for- 
mation shock accele rated CR electrons (Roettiger et "ai][T999at 
lMiniatietalJl2001al) . 

Radio halos and relics are preferentially found in clusters 
showing signs of merger activit\FT. Radio halos are a larger 
mystery than relics. Since the radiating CRe are short lived (0.1 
Gyr, see Fig. |2]i they have to be replenished or re-accelerated in- 
situ. Two classes of models are currently under discussion, the 



e.g. A2256: iBridle & Fomalontl dl976l): iMasson & MaveH fl978); 
Bridle et all <fl97 9): Rottg ering etal] dl994h : Tciarke & ErJflml (2006); 



Kale & Dwarakanath (2010) 



e.g. A 3667 and A3376: ISchilizzi & McAdarrJ dl975l): iGoss et alj 
< 19821); Hemes &McAdaml d!992l) ; iRottgering etal] dl997l) : 
Bagc hietalJd2006l) 

9 Radio gischt in the taxonomy of IKempner et all (2004): 
EnBlin etal] d 19981): iRoettiger et all dl999ah : iKassim et all fcoOll): 
Miniati et all d2001al) : iBagchi et al l d2002h; iBerrington & Dermd 
20031): iHoeftetalJ d2004l 120081): ISiemieniec-Ozieblol d2004l): 
Keshet et alj d2004l) ; iHoeft & Briiggen d2007l) ; iGiacintucci et alj 



of 



Kempneretall d2004h: 
d2001l): lEnfilin etal] 



2008): lBattaglia e t al. (2009) 
10 Radio phoenix in t he tax on omy 
EnBlin & Gopal -Krishnl d2001l); ISleeet al . 

d200ll): lEnfilin & Briiggenl d2002h: iKaiser & Cottej d2002h : 

lBruggenetal.ld2003l);lGardM& Ricked d2004l) 

" iBuotd d200ll): ISchuecker et all d200lh: iMarkeyitch & Vikhlininj 
( 2001 1 ) ; iGovonietal] d2004h : IVenturi et al] d2008l) : ICassano et al] 
( 20ld) 



hadronic and the re-acceleration model. A combination or co- 
existence of their underlying processes might be in operation in 
clusters, as the re-acce lerated CRe could be i njected by hadronic 
CRp-gas interactions dBrunetti & Blasi 2005). However, this re- 
quires some fine-tuning of parameter£3 an d it is more likely that 
one process is responsible for most of the radio halo emission 
and the other process is subdominant. Alternatively, different re- 
gions of the sa me halo could be gene rated by the two differ- 
ent mechanisms. iPfrommer et alj d2008l) proposes that the cen- 
tral part of cluster radio halos is hadronic due to the high target 
density there, whereas at the outskirts shock waves have higher 
Mach numbers and can provide Fermi I acceleration. 

Cluster radio halos come in two sizes: cluster wide and there- 
fore giant radio halos and radio mini-halos. The former are pre- 
dominantly found in clusters showing merger activities whereas 
the latter are found in very relaxed clusters which developed a 
cool core, which harbors the mini-halo. The radio luminosity of 
giant halos seems to be strongly correlated with the X-ray emis- 
sivity of the cluster dLiang et alJl2000t iBrunetti et al.ll2009l and 
see Fig. [TJ. The radio luminosities of the mini-halos also seem 
to correlate in the same way with the cluster X-ray luminosity, 
which itself is usually dominated by the cool core emission. 

A large fraction of clusters do not exhibit significant radio 
halo emission of any kind, and only upper limits to their syn- 
chrotron flux are known. About half of the radio deficient clus- 
ters, for which we have Chandra data, show clear evidence for 
some level of cool core structure (Kq < 50 keV cm 2 ) as can be 
seen in Fig. Q] This could either imply that these clusters are 
in the intermediate state between having giant radio halos be- 
cause of merging activity and having mini halos due to strongly 
developed cool cores. On the other hand there could be two pop- 
ulations of clusters - cool cores and non-cool cores - and the 
corresponding radio luminosity responds sensitively to the level 
of injected turbulence by either AGN or cluster mergers, respec- 
tively. 

Therefore, clusters at the same X-ray luminosity seem to be 
bimodal with respect to their radio luminosity. Either they have a 
prominent halo or they do not exhibit any detectable diffuse radio 
halo emission. This indicates the existence of pronounced and 
rapidly operating switch-on/switch-off mechanisms, which are 
able to change the ra dio luminosity by at least a factor of 10 - 30 
(IBrunetti et al.ll2.009t) . Thus a mechanism is easily realised in the 
re-acceleration model of halo formation, due to the short cooling 
time of the radio emitting electrons, which just cool away once 
turbulence is unable to maintain them. It is less easily realised 
in the hadroni c model via magnetic fiel d decay after turbulenc e 
as proposed bv lKushnir et al.ld2009l) and lKeshet & Lo eb (2010), 
since the turbulent decay take s about a Gyr and the m agnetic 
field decay is relatively gentle (Subramanian et al. 2006). In ad- 
dition, this argument appears weak because /iG strong magnetic 
fields are commonly observ ed in clusters without diffuse radio 
emission (IClarke et alj|2001l) . 

However, the assumption that clusters evolve only vertically 
on the L v - Lx plane, as it is the basis of the above argumenta- 
tion on bimodality, is probably incorrect. The central entropy of 



12 The argumentation goes as following: If secondary CRe from 
hadronically interacting CRp get an energy boost by some factor by re- 
acceleration in order to produce the observed radio halo, a similar boost 
(or even larger) can be expected for the CRp. However, such a boost of 
the CRp population would increase the amount of secondary CRe to the 
level required to explain the radio halo directly. Only if some plasma 
physical reasons can be found why CRe are more efficiently acceler- 
ated than CRp (despite the higher losses of the former), such a hybrid 
scenario can be expected to operate. 
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redshift 






#0 




^1.4GHz 


Reference 






10 44 er 


g s 1 


keV cm 2 


10 31 erg s- 1 Hz -1 




giant radio halos 






1E50657-558 


0.2994 


23.03 


t 1.81 


299.4 


± 19.6 


28.21 ± 1.97 


4, 24 


A209 


0.2060 


6.29 ± 


0.65 


100.7 


± 26.3 


1.19 ± 0.26 


4,28 


A520 


0.2010 


8.83 ± 


1.99 


325.5 


±29.2 


3.91 ± 0.39 


12, 2 


A521 


0.2475 


8.18 ± 


1.36 


201.6 


± 36.1 


1.16 ± 0.1 1 


4,9 


A545 


0.1530 


5.66 ± 


0.49 






1.48 ± 0.06 


4, 2 


A665 


0.1816 


9.84 ± 


1.54 


134.6 


± 23.5 


3.98 ± 0.39 


12, 19 


A754 


0.0535 


4.31 ± 


0.33 


70.4 + 


23.8 


1.08 ±0.06 


11, 2 


A773 


0.2170 


8.10 ± 


1.35 


244.3 


± 31.7 


1.73 ± 0.17 


12, 21 


A 1300 


0.3075 


13.97 


t2.05 






6.09 ± 0.61 


4, 15 


A 1656 (Coma) 


0.0231 


3.77 ± 


0.10 


154.0 


±43.0 


0.72 ± 0.07 


11, 23, 10, 25 


A1914 


0.1712 


10.71 


t 1.02 


63.3 + 


22.3 


5.21 ± 0.24 


11, 2 


A2163 


0.2030 


23.17 


t 1.48 


437.3 


± 82.7 


18.44 ± 0.24 


4, 16 


A2219 


0.2281 


12.73 


t 1.37 


411.6 


± 43.2 


12.23 ± 0.59 


12, 2 


A2254 


0.1780 


4.32 ± 


0.92 






2.94 ± 0.29 


11, 2 


A2255 


0.0808 


2.65 ± 


0.12 


529.1 


± 28.2 


0.89 ± 0.05 


11, 22 


A2256 


0.0581 


3.81 ± 


0.17 


349.6 


± 11.6 


0.68 ± 0.12 


11,8,6 


A2319 


0.0559 


7.40 ± 


0.40 


270.2 


±4.8 


1.12 ± 0.1 1 


11, 15 


A2744 


0.3066 


12.92 


t2.41 


295.1 


± 113.4 


17.16 ± 1.71 


4,21 


CL0016+16 


0.5545 


18.83 


1 1.88 






6.74 ± 0.67 


27, 19 


MACSJ0717 


0.5548 


24.60 


t0.3 


158.7 


± 111.6 


50.00 ± 10.00 


14,31,5 


RXCJ2003.5-2323 


0.3171 


9.25 ± 


1.53 






12.30 ±0.71 


4, 18 


radio mini-halos 






A426 (Perseus) 


0.018 


8.31 




19.4 + 


0.2 


4.40 


1,26 


A2142 


0.089 


10.89 




58.5 + 


2.7 


0.66 


30, 19 


A2390 


0.2329 


13.43 


t3.16 


14.7 + 


7.0 


9.77 ± 0.45 


12, 2 


A2626 


0.0604 


1.96 




23.2 + 


2.9 


0.43 


30, 20 


PKS0745-191 


0.1028 


14.06 




11.9 + 


0.7 


27.00 


30,3 


RXCJ1314.4-2515 


0.2439 


10.94 


t 1.81 






0.75 ± 0.15 


4, 28, 17 


Z7160 


0.2578 


8.41 ± 


2.12 


18.8 + 


3.2 


2.19 ±0.26 


12, 7 


no radio halo detection 






A141 


0.2300 


5.76 ± 


0.90 


144.1 


±31.3 


<0.36 


4, 29 


A611 


0.2880 


8.86 ± 


2.53 


124.9 


± 18.6 


< 0.40 


13,29 


A781 


0.2984 


11.29 


t2.82 






<0.36 


12, 29 


A 1423 


0.2130 


6.19 ± 


1.34 


58.8 + 


12.6 


< 0.41 


12, 29 


A2537 


0.2966 


10.17 


t 1.45 


106.7 


± 19.6 


< 0.50 


4, 29 


A2631 


0.2779 


7.57 ± 


1.50 


308.8 


±37.4 


<0.39 


4, 29 


A2667 


0.2264 


13.65 


t 1.38 


12.3 + 


4.0 


<0.42 


4,29 


A2697 


0.2320 


6.88 ± 


0.85 






<0.40 


4, 29 


A3088 


0.2537 


6.95 ± 


1.20 


32.7 + 


9.5 


<0.42 


4, 29 


RXCJ1 115.8+0129 


0.3499 


13.58 


t2.99 


14.1 + 


5.1 


<0.45 


4, 29 


RXCJ15 12.2-2254 


0.3152 


0.19 ± 


1.76 






<0.63 


4, 29 


RXJ0027.6+2616 


0.3649 


12.29 


t3.88 






<0.68 


13,29 


RXJ1532.9+3021 


0.3450 


16.49 


t4.50 


14.3 ± 


1.9 


<0.62 


12, 29 


RXJ2228.6+2037 


0.4177 


19.44 


fc5.55 






<0.91 


13,29 


S780 


0.2357 


15.53 


t2.80 






<0.36 


4,29 


Z2089 


0.2347 


6.79 ± 


1.76 






< 0.27 


12,29 


Z2701 


0.2140 


6.59 + 


1.15 


34.0 + 


4.2 


<0.42 


12, 29 


Z5699 


0.3063 


8.96 ± 


2.24 






<0.54 


13,29 


Z5768 


0.2660 


7.47 + 


1.66 






<0.36 


13,29 


Z7215 


0.2897 


7.34 ± 


1.91 






<0.55 


13,29 



Table 1: Cluste r sample with rad io halo detections and upper limi ts. Sample base from IBrunetti et alj d2009t) . Four mini-halos 
are added from lGitti et alj d2004l). The X-ray luminosities are as in IBrunetti et al.l d2009l) . for the four additional mini-halos data 
was added from iReiprich & B ohringerl (120021) . and for A2626 from IStott et al] (120081) . Central values for the entropy indicator 
Ko = kT^ n n~% 3 are taken from the extrapolation method in ICavagnolo et al.l (l2009t) applied to Chandra data. Kq of Coma at 12 
kpc is from lRaffertvet all d2008l) . Referenc es: 1 = lAllen et all d 19921). 2 = Bacchi et al. d2003l). 3 = lBaum & O'Deal d 19911). 4 = 
Bohringer et al. d2004l)." 5 = Bonafede et al d2009l) . 6 = Brentjens 72008l), 7 = Cassano et al. d2008l). 8 = Clarke & EnBlin d200r3). 9 
= Dallacasa et al. (2009), 10 = Deiss et a l. d 19971) . 1 1 = Ebeling et al. d 19961) . 12 = Ebeling e t al. dl998l) . 13 = Ebeling et al. d2000|) . 
14 = Ebeling et al d2007l) . 15 = F eretti (l2002h. 16 = Fere tti et al. (l200lh . 17 = Giacintucci d2007l) 18 = G iacintu cci et al. Ho09), 
19 _ G iovannini & Feretti 12000), 20 = iGitti et alj d2004l). 2 1 = Govoni et al d2001bl) . 22 = Govoni et al. d2005l) . 23 = Kim et al. 
dl990|) . 24 = Liang et al. d2000h 25 = Rafferty et al. d2q08l) . 26 = Sijbring d2007h.27 = Tsuru et al. (1996), 28 = Venturi et al. 
120071) . 29 = I Venturi et al.l d2008l) . 30 = I White et all d 19971) . 31 = van Weeren et al. d2009l) . 
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Fig. 1: Correlation of radio halo luminosities with cluster properties of clusters in Tab. Q] Left: Radio halo luminosity vs X-ray 
luminosity. Right: Radio halo luminosity vs central entropy indicator Ko for the subsample of clusters for which high resolution 
Chandra data are available. 



many of the radio halo deficient clusters is low, indicating that 
they host or are forming a cool core. Cool cores tend to dom- 
inate the X-ray luminosity of clusters. Therefore, also a strong 
horizontal evolution in this plane can be expected once the radio 
halo luminosity decreases as the cluster relaxes after a merger. 
Radiative cooling in the central cluster regions increases the 
gas density and therefore the X-ray luminosity while it also de- 
creases Ko. This alleviates the requirements on the speed and the 
magnitude a radio halo switch mechanism must fulfill in order 
to explain the observations. 

Nevertheless, such a mechanism is probably needed to un- 
derstand the observational data. It has to full-fill two require- 
ments: it should be able to extinguish a radio halo, but also 
not prevent it from being switched on again later, otherwise we 
would not observe halos in the present universe. The mechanism 
we are discussing in this work is actually able to switch radio ha- 
los on and off in both models, the hadronic and re-acceleration 
model. Since both have their individual strengths and weak- 
nesses in explaining the different observational features of radio 
halos, we want to briefly discuss those first. 

2.2. Hadronic models 

In the hadronic model the accumulated CRp inject continuously 
the radio emitting CRe into the ICM due to well known hadronic 
process pcr + p — > 7T* + . . . — > e ± + v e /v e + + + . . . , 
(see footnote Q] for references). The CRe loose their energy 
nearly instantaneously (see Fig. [2]) and therefore the radio emis- 
sion traces a combination of the CRp population, the distribu- 
tion of the target gas density, and the magnetic field profile of 
the galaxy cluster in the case of weak magnetic field strengths 
B < Bcmb = 3.27 (1 + z) 2 pG (z being the redshift) and without 
magnetic field dependence in the strong-field case. 
The hadronic model has a number of advantages: 

- Since CRp are expected to be present in the ICM due to 
structure formation shock waves, active galactic nuclei, and 
the deposition of interstellar media of galaxies, radio halos 
emerge very naturally in these models. 



- Also the observed power- law spectra of ma ny radio halos 
(e.g. in the Bullet cluster, Liang et al. 2000) emerge natu- 
rally, since they reflect the power-law spectra of CRp accel- 
eration. 

- The observed smooth, regular morphology of halos is also 
expected in hadronic models, since the long lived CRp had 
enough time to become distributed within the cluster volum e 
(Mini ati et al.ll2001bHMiniatill2003l:lPfrommer et al.ll2008l) . 

- The observed correlation of radio halo luminosity with X-ray 
luminosity (L v -Lx relation, see Fig. UJ of clusters with halos 
holds within these models for plausibl e magnetic field val- 
ues and CRp acceleration efficiencies (Dolag & EnBlin 2000; 
iMiniati et al.ll2001al:|Pfrommerll200llKeshet & Loebll 2010). 

- Also the correlation of radio halo surface brightness with X- 
ray surface brightness of clusters with halos holds roughly 
for sensible magnetic fi eld profiles a nd CR p acceleration 
efficiencies. [However, iGovoni et al.l (1200 lal) find that the 
naively expected radio halo profile does not seem to fit the 
observed asymptotic in case of weak magnetic fields (B < 
3/*G).] 

However, there are also a number of issues with the hadronic 
model: 

- About two thirds of the most X-ray luminous clusters do not 
exhibit radio halos, whereas the hadr onic model seems to 
suggests that all clusters exhibi t halos (|Miniati et aLlu OOla; 
Miniatil 120031 |Pfrommerll2008t iPinzke & Pfrommerfe oiO). 
[This problem will be addressed and alleviated by this work.] 

- The curvature in t he total spectrum cl aimed for the Coma 
cluster radio halo dBrunetti et al.1 12001) are not reproduced 
in current numerical models of the hadronic scenario. 
[However, the particle transport in these models neglects dif- 
fusion and streaming which is potentially important. In addi- 
tion, the high frequency observations of the steep spectrum 
radio halo, on which this claim rests, are extremely diffi- 
cult, since point sources much brighter than the diffuse radio 
halo itself have to be subtracted, and the sensitivity to large 
scale emission of radio interferometers is reduced at higher 
frequencies. Thus the curvature might be partly a fluke, as 
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also some inconsistencies of flux measurements at the same 
frequencies in thi s spectrum indicate (s ee 0.4 and 1.4 GHz 
fluxes reported in Thie rbach et al.ll2003l) .l 

- Spectral ste epening at the edges as proposed by s ome ob- 
servations (Br unetti et alj|2001t iFeretti et al.ll2004al) can not 
be explained. [However, the observational arguments above 
are even more severe for this. In addition, the negative 
flux of the Sunyaev-Zel'dovich decrement of clusters af- 
fects the spectrum especially hard in the outs kirts, and also 
provides some bending to the total spectrum (EnBlin 2002; 
iPfrommer & Enfilinll2004bl) . Finally, energy dependent CRp 
transport would create such spectral variations, as argued in 
this work.] 

The hadronic model makes one hard prediction, which hope- 
fully will permit its confirmation or rebuttal at some point in 
the future. The radio halo emission should always be accom- 
panied by some level of gamma-ray flux, due to the hadronic 
production of neutral pions and their decay into gamma-rays, 
Pes. + P — > 7i° + . . . — > 2y+ ... The current upper limits on dif- 
fuse gamm a-ray flux from cluster of galaxies by the Fermi col- 
laboration (lAckermann et alj|2010 ) are still well above the pre- 
dictions of expected fluxes, even for the mos t optimistic assump- 
tions about the CR acceleration efficiency (Pinz ke & Pfromrneil 
1201 Ol) or by tying the expected gamma-ray emission to the simu- 
lated radio halo emission for rea sonable assumed magnetic field 
strengths (Do nnert et al.l l2010aUbl) . They are far off the mini- 
mal gamma-ray flux ex pected in the lim i t of strong magneti c 
field strength ( » 3 uG: iPfrommerl l2008t lAleksic et alJl2010h . 
iBrunettil d2009l) and iJeltema & Profumol OOlOl) argue that un- 
usually strong ICM magnetic fields would be required by the 
hadronic model for radio halos with reported steep radio spec- 
tra (e.g. Abell 1914, 2256). However, this is only true if the 
CRp spectra can be extrapolated from the radio emitting energies 
(~100 GeV) into the subrelativistic regime (< 1 GeV) without 
any spectral break as assumed by these authors. 

2.3. Re-acceleration models 

In re-acceleration models, a pre-existing CRe population at 
lower energies of about 0.1-10 GeV gets re-accelerated into 
the radio emitting regime of about 10 GeV by plasma waves 
(see references in footnote [2J. These are generated by the tur- 
bulence during and after a cluster merger event. Some level of 
re-acceleration has to happen most of the time or frequently 
enough in order to prevent the CRe population in the cluster cen- 
ter from loosing its energy completely due to Coulomb losses on 
a timescale of about 1 Gyr. 

Also the re-acceleration model has its advantages: 

- The bimodality of radio halo luminosities is explained in 
this model by the presence and decay of the re-accelerating 
turbulence in merg ing and relaxed clusters, respectively 
dBrunetti et al.ll2009l) . 

- CRe are expected to be accumula ted in the ICM due to 
injection by radio galaxy outflows dGiovannini et al l l 19931) 
and acceleration at shocks, if some level of continuous re- 
acceleration can prevent them from thermalisation. 

- Diffusive re-acceleration is a natural plasma process, which 
must occur in turbulent astro physical environments as clus- 
ters (ISchlickeiser etalJll987h . 

- The complex morphologies reported for some radio halos 
come natur ally about due to the effect of intermittency of 
turbulence dFeretti et al.ll2004al) . 



- Significantly curved radio spectra are also very natural, since 
the interplay of ac celeration and cooling produces spectral 
bumps and c ut-offs dFeretti et al.l2004at IBrunetti & Lazarianl 
120071 l201Ch . 

- The Ly—Lx relation h olds also for reasonable magnetic fields 
dCassano et alj|2007l) . 

The issues with the re-acceleration model are: 

- Second order Fermi acceleration is known to be very ineffi- 
cient since the efficiency scales with (u W ave/c) 2 <k 1. The 
acceleration efficiencies used in the re-acceleration mod- 
els are difficult to be derived from first prin ciples and are 
often fit to reproduce the observations (e.g. iBrunetti et al.l 
(2001) adopt the magnetic field profile to match their re- 
acceleration profile). [However, magnetosonic turbulence 
may be more efficient, and could provide sufficient re- 
accelerati on if the compressible turbulence is of the required 
strength. IBrunetti & Laza rian ( 20071 120101) show that a com- 
pressible wavefield of 15-30% of the thermal energy con- 
tent of the cluster, reaching down to small scales with a 
Kraichnan-scaling, would be sufficient. It should be noted 
here that Fermi I acceleration at ICM shock waves is also a 
promising re-acceleration mechanism.] 

- The maximal cooling time of CRe is r e ~ 1 Gyr(« e /(3 x 
10~ 3 cm 4 )) (see Fig.[2j». Without replenishing the 100 MeV- 
seed population of CRe, the re-acceleration would not be 
able to generate 10 GeV CRe that radiate GHz radio waves. 
Either continuous injection or re-acceleration of the short- 
lived central CRe population is needed to counteract cool- 
ing or a process that moves CRe inwards without much 
losses from the dilute outer cluster regions where the cool- 
ing time approaches the Hubble time. [This problem will be 
addressed and alleviated by this work.] 

- Any re-acceleration mechanism is unable to discriminate be- 
tween CRe and CRp at GeV energies. In case re-acceleration 
operates, the CRp population gains more than the CRe pop- 
ulation in the long run, due to the much lower losses. Given 
that the injection efficiency of CRp by shock waves is ex- 
pected to be hundred times higher than the CRe one, the re- 
accelerated CRe can easily be outnumbered by hadronically 
injected CRe. [However, the possibility exists that the CR 
population is dominated by injection of radio plasma, which 
might contain mostly CRe.] 

- The natural curvature of re-acceleration radio spectra re- 
quires fine tuning to reproduce the observed power-law radio 
spectra seen in some radio halos. [However, a detailed statis- 
tics of the radio halo spectral slopes that includes statistical 
and systematic uncertainties is not published yet.] 

- The regular morphology of many radio halos might also be 
a challenge to the re-acceleration model, since turbulence 
and therefore re-acceleration is expected to be intermittent. 
[However, detailed numerical simulations of the expected 
morphologies are not yet published]. 

The hardest prediction of the re-acceleration model so far 
is that low X-ray luminous clus ters should not exhibit radio 
halos dCassano et alj|200a, 12008). This prediction will become 
testable with upcoming sensitive radio telescope arrays. 

3. Cosmic ray transport in galaxy clusters 

3. 1 . Cluster weather conditions 

The atmospheric conditions in galaxy clusters differ strongly de- 
pending on the time elapsed since the last cluster merger event. 
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Fig. 2: Cooling time of CR in the ICM as a function of their kinetic energy. Left: CRe cooling times for typical densities and 
magnetic field streng th that range fr om t he central to the peripheral r egions in galaxy clusters. Coulomb and IC/synchrotron cooling 
is modeled following Go uld! dl972bl) an dlRvbicki & Lightmanl (fl979). respectively. Right: CRp cooling times for the same densities. 
Coulomb cooling is modeled following Gouldl (Q~972a). The hadronic cooling time above the kinematic threshold for pion production 
is T pp = 1 / (0.5 (Tpp n N vcr) with cr pp = 32 mbarn, an inelasticity of ~ 50% and uqr the CRp velocity. It is apparent that CRp above 
10 GeV have livetimes in the ICM at least 60 times longer than CRe at any energy. CRe can survive for a Hubble time without 
re-acceleration only within the dilute outskirts of clusters. The radio emitting electrons have an energy of about 10 GeV in fiG 
fields, and herefore a lifetime of 0.1 Gyr or less. If they are of hadronic origin, their parent CRp had energies of about 100 GeV, 
which have considerable longer lifetimes. 



During a cluster merger, bulk velocities close and above the 
sound speed of c s as 1000-2000 km/s are injected on scales of 
a few xlOO kpc. A fair fraction of this energy goes into tur- 
bulent motions, whi ch can account for 10-20% of the the rmal 
pressure of the ICM dSchuecker et al.ll2004t fPaul et alj|2010h . A 
good fraction, however, is directly dissipated in shock waves, 
mostly in the cluster center, however with the highest Mach 
number shock waves and therefore the most efficient CR acceler- 
ation sites appearing in the c luster outskirts (Miniati et al. 2000; 
iPfrommer etani2006L[2008h . 

The cluster turbulence persists for about a Gyr, partly be- 
cause the turbulent decay of eddies takes several eddy turn over 
times to transport the kinetic energy to dissipative scales, partly 
because the gravitational drag of the merging dark matter halos 
continues to stir turbulence for some timq3 

As the turbulence decays, the cluster atmosphere settles into 
a stratification with the lowest entropy gas at the bottom and 
higher entropy gas as larger radii. The turbulence may become 
more two-dimensional, since all radial motions are working 
against gravitational and/or pressure forces. If the cluster be- 
comes sufficiently quiet, thermal instabilities might set in due to 
heat flux along magnetic field lines, rearranging the field in pref- 
erentially radial orientation as predicted by the magneto-thermal 
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instability (IBa lbus 2000; Parrish & Stone] 12007b iParrish et al.1 
120081; IPfrommer & Dursill2010l) . 

At the dense cluster center, X-ray cooling could become 
catastrophic due to the onset of a cooling instability. However, 
the developing cool cores seem to be stabilised against a com- 
plete collapse by energetic f eedback from the ce ntral galaxy in 
form of AGN outflows (e.g. Churazov et al. 2002) or supernovae 
driven winds, both fed by condensing gas from the cooling re- 
gion. Since cool cores have temperatures of only a few keV, their 
sound speed is below 1000 km/s. The turbulence, which is most 
likely subsonic there, will therefore only have a speed o f a few 
100 km/s there dFabian et al.ll2003HEnfilin & VogjboOfih . 

Let us model the cluster weather in a simplistic fashion, in 
order to estimate the order of magnitude of the transport pro- 
cesses. We assume that the turbulence is injected on a length 
scale L tu comparable to the cluster core radius r c ~ 200 kpc, 
which is also roughly the atmospheric scale height in the inner 
part of the cluster: 



L t u =Xtur c with %tu ~ 1- 



(1) 



Although the turbulence injection is supersonic, shock waves 
rapidly dissipate kinetic energy until turbulent velocities u tu be- 
come subsonic. Thus, the turbulence stays nearly trans-sonic for 
about a Gyr. We therefore write 



v tu = am c s with a tu (f) < 1 



(2) 



clearly being time dependent. Note that for notational consis- 
tency reasons, we identify the turbulent Mach number by M tu = 
Ctu- Half an eddy turn over time 



2 v tu 2 a tu 



(3) 
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is probably larger, but comparable to the sound crossing time of 
an eddy t s = L tu /c s . 

The hydrodynamical turbulence is expected to drive a small- 
scale magnetic dynamo, which is believed to saturate for mag- 
netic energies, which are a fraction a\ ~ 0.03 - 1 of the ki- 
netic energy density, and a magnetic coherence length As which 
is smaller than the turbulent injection scale by some f actor 
0(1 - 10) (ISubramanian et al J 12006: EnGlin&Vogj[2006). We 
have to stress that for filamentary fields, the perpendicular scale 
of the magnetic structures, which dominates this number, is 
much smaller than the parallel scale. Flux tubes can have ex- 
tension Lg comparable to the eddy size. 

We parametrise the Alfvenic velocity by 

v A = a B u to = a B arm c s (4) 

and the magnetic bending-scale by 

L B =Xi}L m with xb ~ 1. (5) 

The ICM magnetic field strength is a relatively poorly known 
number, despite substantial efforts to measure it in the different 
cluster environments. Inverse Compton based methods[3 suggest 
rather low field strength of <9(0.3 fiG). These methods identify 
reported excess emission above the thermal ICM X-rays in the 
extreme ultraviolet and hard X-ray bands with CMB photons 
up-scattered by the radio halo emitting electrons. However, the 
identification is not unique, and even the detec tion of the ex- 
cess emission is heavily debated (Molendi 2008, and references 
therein). Recent work u sing the Swift/Burst Alert Telescope 
(Aiell o etaDl2009l l2010l) was able to reject the hypothesis of 
non-thermal power-law emission in a sample of 20 clusters, ex- 
cept for the Bullet cluster where they seem to require a non- 
thermal power-law component. They are able to obtain excellent 
fits of their data with multi-temperature models - hence inferred 
values of the magnetic field strength using the inverse Compton 
method are better regarded as lower limits on the field strength. 
Magnetic field me asurements using Faraday rotation data of 
galaxy clusters (e.g. lClarke et al.l20 01) find fields 0(3 fiG), how- 
ever, these numbers depend on the assumed characteristic field 
scale As- Note that the IC-based magnetic field strengths are 
weighted by the CRe density at ~ 10 GeV while the Faraday 
rotation-based values are weighted by the thermal electron den- 
sity. In the case of the Hydra A cool core, where high quality 
data enables measuring the magnetic power spectrum in detail, 
magnetic fields of (9(10 - 30yuG) see m to support the idea of a 
saturated field strength with a\ ~ 0(1) iKuc har & EnBlin 2009). 

Anyhow, it will turn out that the precise magnetic field 
strength is less important for our argumentation than the char- 
acteristic bending length and the assumption that the magnetic 
energy density is proportional to the turbulent one, irrespective 
the numerical factor in between them. 

14 During the kinematic stage of the small-scale or fluctuation dy- 
namo, the perpendicular scale is the microscopic resistive scale. How 
this changes due to nonlinear effects is n ot yet fully clear (For pa - 
pers on this aspect seelHaugen et alj J2004): Schekochi hin et alJ f2004): 
Brandenburg & Subramanian 12005)) . An increase i n the effective re- 
sistivity due to nonlin ear effects ([ Subramanian 1999) or due to micro- 
scopic plasma effects (Schekochihin & Cowley 2006) is needed if the 
perpendicular scales have to become a fraction of the turbulent scale, 
as required to explain the observed coherence of cluster magnetic fields 
(EnBlin & Vogt 2 00g). 

15 e.g. fRephaeli et alj dl994l); lEnfllin & Biermannl d!998l>: 
Saraz in & LieuNl998h : iBowver & Berghoferl dl998[> : | V61k & Atovanl 
dl999[) : ISarazir] d!999h : ISarazin & KempneJ 12000) : Utovan & VoM 
(2000) 



3.2. CR transport 

The discussion in this section is applicable for CR electrons as 
well as protons. We describe our CR population by their phase 
space density f(r,p,t) = d A N l(dr dp) assuming that the mo- 
mentum distribution is nearly isotropic and can hence be de- 
scribed by a one-dimensional spectrum. Anisotropies lead to 
CR streaming with a velocity v cr ~ c s , as we will argue below. 
Therefore anisotropies are expected only on the order u cr /c CI ~ 
c s /c ~ 10~ 3 which we can safely neglect here except for their 
transport effect. Here c cl < c denotes the CR particle velocity, 
which is close to the speed of light c. We express the CR mo- 
mentum p in units of mc, with m being the particle mass. For 
CRp, we assume a power-law particle spectrum, 

f(r,p,t) = C(r,t)p- a . (6) 

Here a(r,t) m a » 2.1 — 2.5 is the dominating CRp spectral 
index around 100 Ge V expected in the ICM from structure for- 
matio n shock waves dMiniati et alJl200Tbt iPinzke & Pfrommer 
1201 Ol) and C(r, t) the spectral normalisation constant, which will 
capture the spatial variations and temporal evolution. 
The full CR transport equation reads 

d l* (Vin d (pn 3(3 f ) +q -±. (7) 
at an dp an \ orj ) ri oss 

Here, v - v tu + v st is the CR transport velocity, which is the 
sum of the (turbulent) ICM gas velocity, v lu , and the CR stream- 
ing velocity with respect to the gas u st , to be discussed below. 
For relativistic protons, the continuous momentum loss or gain 
term p(r, p, t) can be approximated to be solely due to adiabatic 
energy changes, since Coulomb and radiative losses can be ig- 
nored and hadronic losses are modelled as catastrophic losses 
with a timescale ti oss ~ l/(o- vp n p c) ~ <9(10Gyr) (see Fig. [2j. 
Thus 

p*p° d = ~£v-v. (8) 

The adiabatic energy changes of CRe follow the same equation, 
however, radiative losses and Coulomb losses can usually not be 
ignored for CRe (see Fig. Since we do not want to enter a 
detailed modeling of the CRe spectrum for the sake of briefness, 
we ignore the cooling in our calculations and can therefore only 
make timescale arguments just to see if and how CRe transport 
processes might be relevant. 

The total transport velocity v — v tu + v s( appears in the 
adiabatic term, since v represents the combined frame of the 
plasma waves and magnetic field line motions which together 
confine the particles a bit like walls. Thus, the energy transfer 
from streaming CR into waves is an adiabatic process for the 
CR, the particles push these walls away while expanding into 
unoccupied regions. The energy deposited in the waves gets dis- 
sipated by Landau and other damping processes and the inter- 
play of excitation by the CR streaming and damping determines 
the precise value of the streaming velocity v st . This will be dis- 
cussed more in detail in Sect. 13.31 For the time being we antici- 
pate the result adopted here that 

Vgt = ffstc s with ar s , < 1. (9) 

In addition, we will assume that the diffusive transport (dis- 
cussed below) is possibly less efficient, since if it leads to sig- 
nificant CR bulk motions, these will excite plasma waves which 
limits the transport to the streaming regime. 

The CR diffusion tensor is highly anisotropic with k\\ and k ± 
being the diffusion coefficients parallel and perpendicular to the 
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local magnetic field direction, respectively. Nondiagonal terms 
that mix the parallel and perpendicular direction appear in case 
of diverging or converging magnetic flux. Very often, k\\ » k ± 
so that we can ignore cross field diffusion on short timescales. 
Although several theoretical workfl claim that in turbulent en- 
vironments k ± might become comparable to k\\, the observation 
of sharply edged boundaries of radio plasma questions if the nec- 
essary conditions for this are realised in the ICM. 

When the CR are tied to their flux tube it makes sense to 
introduce an affine coordinate x along the tube and a CR density 
per infinitesimal magnetic flux d 2 4> = Bdydz with g(x, p, t) = 
d 4 N/{d 2 (p dx dp). The transport equation along the flux tube can 
be written as 

dCesarskv& Volklll978b lEnfilinl 12003k . where we dropped the 
injection term q and the catastrophic losses. The magnetic field 
enters the transport equation since magnetic mirrors reflect par- 
ticles away from locations with converging field lines, so that 
the CR population tries to establish a constant density per vol- 
ume, and not per magnetic flux. Regions where the magnetic 
field strength is lower and therefore the flux tube diameter is 
larger will carry more CRs per tube length. Since for a flux tube 
which connects volume elements on the bottom and top of an 
atmospheric scale height (or the fraction ;^;^tu of it) the tube di- 
ameter is expected to be larger at larger cluster radii. This is due 
to the low-pressure environment which enables the gas to expand 
there and CR will tend to occupy preferentially more peripheral 
regions. 

3.3. CR streaming velocity 

Cosmic ray streaming is an important ingredient of the present 
work, as mentioned above. We discuss here briefly the range of 
possibilities for the CR streaming speed. But before we enter the 
detailed discussion, we summarise our main insights: 

For a very turbulent medium, CR may be well trapped by 
magnetic irregularities produced by the turbulence. For a less 
turbulent high-/? plasma, the sound speed seems to be a ref- 
erence speed However, in case the turbulence and CR en- 
ergy densities are low, even higher transport speeds are possi- 
ble. Since we only know the streaming speed approximately, 
we will parametrise it in terms of the sound speed. Note, that 
complex magnetic topologies reduce the effective, macroscopic 
speed, and this effect should be stronger for turbulent clusters 
than for relaxed ones. 

Now we investigate the different factors determining the 
streaming speed. In doing so we shall clearly distinguish be- 
tween two separate regimes, one for merging clusters and the 

16 [jo kipii ( 196 &ll967h:lRechester & Rosenblutrj dl978h:lDuffv et all 
dl995h: lBiebeT& Matthae usril997l): iMichalek & Ostrowskil dl997L 
1998): lGiacalone & Jokipiil dl999T )- |Naravan & Medvedevl d200lh 

17 In a low-/? plasma, the CR streaming velocity is linked to the Alfven 
velocity, which exceeds the sound speed there. However, this can obvi- 
ously not be true in a high-/? plasma. It would imply that in the limit 
of vanishing magnetic field strength the CR get completely immobile 
due to the vanishing Alfven speed. Already intuition tells us that this 
must be wrong, since for disappearing magnetic fields, the coupling of 
CRs to the plasma gets weaker and therefore the CRs should stream 
faster. Thus, there must be a characteristic velocity, below which the 
Alfven velocity is not limiting the streaming velocity any more. Plasma 
physical arguments indicate that this is roughly the sound speed. 



other for relaxed clusters. In the first case the ICM is highly tur- 
bulent. The turbulence is transonic, i.e. the sonic Mach number 
is M = a tu > 1 and, for a typical magnetic field strength of 
a few yt/G and sound speed of 10 3 km/s, highly super-alfvenic, 
i.e. the Alfvenic Mach number is Ma = v^/vp, ~ 10. A small- 
scale dynamo is also operating at this stage. We assume that 
this has amplified the field to the observed level of a few mi- 
crogauss from perhaps weaker fields. Under these conditions the 
magnetic field lines are easily bent by the inertia of the ICM 
motions, at least on scales where the turbulence remains super- 
alfvenic. If the turbulence is injected on scales L tu a few x 10 2 
kpc, and follows Kolmogorov's scaling, then super-alfvenic con- 
ditions persist down to scales La — L tu M A 3 a few x 10 2 pc. 
Under these conditions the CR particles are efficiently trapped in 
between large amplitude magn etic field fluctuations by adiabatic 
mirroring (e.g. iKulsrudll2005l) . So, at least macroscopically, the 
CR particles effectively follow the advective motion of magnetic 
field lines. 

As the merger process approaches dynamical equilibrium, 
the amount of ICM turbulence drops considerably in the core 
regions, with t ypical values of the gas speed u tu ~ a few x 10 2 
km/s (see, e.g. Faltenb acher et al.ll2005t IPaul et alJ l2010). This 
changes the dynamical picture in the core regions at a qualitative 
level. In fact, the turbulent motions are now subsonic and, more 
importantly, trans-alfvenic, i.e. Ma > 1. Under these conditions 
adiabatic mirroring will be modest and the most effective way to 
limit the propagation speed of the particles is through the self- 
generated turbulence. 

In fact, CRs in the Solar neighborhood are observed to have 
a very small anisotropy (about 1 part in 10 4 ). This is explained 
by pitch-angle scattering of CRs by plasma waves, in partic- 
ular Alfven waves, in our interstellar medium. These waves 
can be generated by the CRs themselves, by a streaming in- 
stability, as they gyrate around and stream alon g the magnetic 
field at a velocity faster than the Alfven speed (|Wentzell ll968: 
iKulsrud & Pearcdll969tlKiilsrud & Cesarskvlll971h . The waves 
traveling in the direction of the streaming CRs (and satisfying a 
resonance condition for low frequency waves, £ wave y cr v cs: fi a - = 
Q cr ) are amplified by this streaming instability. Here £ wave is 
the wavenumber of the amplified wave, v a is the CR velocity, 
y cr its Lorentz factor, O cl its gyro frequency, and p CI = cos cr , 
with Oct the pitch angle, the angle between the CR velocity and 
the magnetic field direction. For p CI ~ 1, the amplified waves 
have a wavelength of order the CR gyro-radius. These self- 
generated waves can then resonantly scatter the CRs and lead 
to an isotropisation of the pitch angle distribution in the wave 
frame. Hence in the low plasma-/? interstellar medium the CRs 
are expected to stream wit h respect to the gas with the wave (or 
Alfven) velocity (see, e.g., Kulsrud 2005, for a detailed discus- 
sion). The scattering also leads to a transfer of momentum and 
energy from CRs to the waves, which is in turn transmitted to the 
thermal gas, when the waves damp out, although this is a higher 
order effect than the pitch-angle scattering process and therefore 
slower. 

A potential problem with the above picture is due to the fact 
that in order to isotropise their pitch angle distribution some CR 
particles need to be scattered by the waves to have negative p cr , 
or reverse their parallel component of the velocity. This can only 
happen if they get scattered through p CI = 0. But for fi cr ~ 0, 
the resonance condition implies that resonant waves will have 
a large £ wave . In a high-/? plasma, such as the ICM, such large 
fc wave waves can suffer str ong resonant damping by thermal pro- 
tons dHolman et al.ll 1 979b - Thus to maintain them at a significan t 
level one needs a larger streaming speed. Hol man et al.l ([l979) 
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then deduced that, in a high-/? plasma, CRs can then stream at a 
speed of order or greater than the ion sound speed, which in this 
case will be larger than the Alfven speed. 

Possible processes can help in crossing t he yu cr - gap, in- 
clude mirroring on lon g wavelength waves dAchterberdl 19811 
iFelice & Kulsrudl [2001), an d resonance broadenin g which re- 
laxes the resonant condition (Achterberg 19 811 lYan & Lazarianl 
l2008h . As for the former case, long wavelength waves - respon- 
sible for the mirroring - are assumed to be generated by the 
streaming CR particles themselves, albeit at a smaller pitch an- 
gle. However, they can also be produced as a result of resid- 
ual turbulent mo tions which we know exi s t even in the IC M of 
relaxed clusters dFaltenbacher et all 120051; iPaulet all 1*20 1 Oh . In 
any case lFelice & Kulsrudl (1200 lh find that the contribution from 
mirroring to the diffusion coefficient is important but not domi- 
nant. The streaming velocity, v str , can then be found by balanc- 
ing (A) the growth rates of the resonant waves produced by the 
streaming particles and the (nonlinear) damping rates, which is 
proportional to the amount of resonant waves; and (B) by noting 
that the streaming velocity is determined by the anisotropic com- 
ponent of the CRs which is inversely proportional to the amount 
of resonant waves responsible for pit ch angle scattering. Thi s 
is basically the argument followed bv lFelice & Kulsrudl d200ll) . 
If we apply this argument to the case to the ICM we find that, 
among others, the streaming velocity scales inversely with the 
square-root of the number density of the CR themselves. If the 
CR energy is at the level of ten percent or so, then the streaming 
velocity is of order of the thermal speed of the ions, i.e. we find 
that Eq. [9] is approximately valid. If the CR population is sub- 
stantially lower than that, e.g. as in the case of CRe or for CRp 
after a result of their streaming motions, the streaming velocity 
can become significantly larger than the sound speed. 

We should note that during both the turbulent and quies- 
cent phases, the magnetic field is not static. When the gas is 
turbulent, there is a dynamo, which not only amplifies the field 
but also makes it spatially and temporally intermittent. A place 
which has strong field at some time will, after eddy turn over 
time, have weaker field and also possibly with different con- 
nectivity. However this does not affect the streaming instability 
drastically as the growth and damping rates of the waves will 
generally occur at a much faster rate than the growth rate of the 
small-scale dynamo. But it will affect the long time behaviour of 
CRs, because of the constantly changing and complex connec- 
tivity of the field, although the consequences of such changes 
for CR propagation have yet to be studied. More interesting is 
what would happen when there is no longer merger-driven turbu- 
lence. The field itself then generates decaying MHD turbulence, 
with the kinetic and magnetic energy decaying as a power law 
in time, and th e coherence length of th e field increasing also as 
a power law dSubramanian et al. 2006). Such an increased cor- 
relation length and simpler connectivity could make it easier for 
the CRs to propagate throughout the cluster, albeit limited by the 
streaming velocity. 

3.4. Interplay of advection and streaming 

First, we consider an isolated magnetic flux tube with CRs con- 
fined to it to illustrate the interplay of advection and streaming 
with a basic picture. This will be generalised later on in Sect. 13.61 
to a more complex situation including CRs escaping from mag- 
netic structures. Focusing on the CR transport along a single flux 
tube embedded in the ICM, we start with the idealised picture of 
a static flux tube frozen into a static plasma as shown in Fig. [3] 
on the left. Any central concentration of CR will escape due to 



CR streaming on a timescale of 

Lb Xb L t u n 
r st = — = . (11) 

This leads to a homogeneous CR distribution within the flux tube 
(Fig. [3] middle). Turbulence turns the magnetic structure upside 
down on half an eddy turnover time r tu = 7rL tu /(2 v tu ). This is 
comparable to, or less than, the CR escape time, 

^yn.-**— (12) 

and thus a good fraction of the CR from larger radii will be com- 
pressed towards the center, from where they again start stream- 
ing to larger radii. 

The transonic turbulence is therefore able to maintain a cen- 
trally enhanced CR density by pumping expanded CR popula- 
tions downwards. As soon as the turbulent velocities become sig- 
nificantly subsonic, this pumping becomes inefficient, since the 
streaming will be faster than the advection. At this point a nearly 
constant volume density of CR establishes within a closed flux 
tube, meaning that most CR are residing at larger cluster radii. 
However, the CR density can again become centrally enhanced, 
if the turbulence becomes strong again, e.g. during the next clus- 
ter merger. 

During the phase of transonic turbulence, we regard two 
Lagrangian volume elements which are exchanged radially by an 
eddy, one starting at small radius with volume V\ and the other 
at large radius with volume V2 — XV\, with X > 1. We assume a 
relaxed CR population with power law spectrum as in Eq. |6]be- 
ing present in them, with C\ — C2 = C. Now, V2 — > V' 2 = V\, the 
enclosed CR get compressed and gain momentum adiabatically 
according to 

P2^>P 2 = X ~P2- (13) 

The compressed CR population is again described by Eq. [6] but 
now with 

C' 2 = CX a -^. (14) 

Correspondingly, the expanded CR population in V[ has a 
power-law normalisation constant of C\ = C x^ a+2) ^. A simple 
swap of the two volumes can generate a CR density contrast (at 
fixed momentum) of C' 2 IC[ = x <2a+4 ^ 3 w 2 3 = 8 from a previ- 
ously homogeneous distribution if the ratio of the initial volumes 
was only X = 2. 

The complete process of advective compression and expan- 
sion via CR streaming is, from the perspective of the CRs, com- 
pletely adiabatic. The energy losses during the streaming ex- 
pansion phase will be exactly compensated during the advec- 
tive compression phase. This might surprise, since the dissipa- 
tion of the excited waves is certainly not a reversible process 
(for any practical purpose). But the energy injected into the 
wave fields comes from adiabatic work done by the expanding 
CR gas. All what the CRs feel are the expanding walls of the 
self-generated plasma waves, and the converging magnetic fields 
confining them when they are dragged back inwards by a down- 
wards directed flow. 

However, from the perspective of the turbulence, the process 
is dissipative. Not all of the kinetic energy invested into CR com- 
pression will be returned to the kinetic energy budget during a 
gas expansion phase, but some fraction is transferred via plasma 
wave generation and damping into thermal heat. Thereby this 
process contributes to the damping of turbulent flows as well as 
to the heating of the ICM. Since the CR can release their en- 
ergy at different locations than where they got it via adiabatic 
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Fig. 3: Sketch of the interplay of CR streaming and turbulent advection for a single flux tube in a stratified atmosphere with gravity 
pointing downwards. Left: The dense CR at the center stream along the tube towards the CR depleted regions at larger atmospheric 
height. Middle: CR streaming stops as soon as a homogeneous CR space density is achieved. A turbulent eddy (represented by 
its angular momentum axis) starts to turn the magnetic structure upside down. Right: The former outer parts of the flux tubes 
are compressed at the center, and harbor now an overdense CR population, whereas the former inner parts are expanded at larger 
atmospheric scale height and therefore have now an underdense CR population. Again CR streaming sets in. 



compression in the turbulent flow, the spatial distribution of their 
heating power will in general differ from that of the turbulence. 
Since CR actively propagate preferentially to regions of low CR 
densities, their heating profile can be expected to be more regular 
than that of the probably intermittent turbulent cascades. 

3.5. Diffusive acceleration 

CR diffusion will always be present at some level, although it 
is probably a less efficient transport mode compared to stream- 
ing. It has, however, some other advantage in that this propaga- 
tion mode does not lead to adiabatic energy losses of the CR. 
Diffusion also tends to establish a constant CR density per vol- 
ume as streaming does, but with a positive energy balance if 
counteracted by turbulent advection. If the volume expansion 
factor X can be split into a (dominant) streaming and (subdomi- 
nant) diffusion part, we obtain 

X=X stx X<xff. (15) 

If this is identical to the (inverse of the) subsequent advective 
compression X^, then a net energetisation of the CR population 
happens. This is described by 

g-l 

^*C^L (16) 

dt T tu 

and leads to an exponential energy gain of the CR. The individ- 
ual particles gain per cycle the momentum 

p^p'=pXl s . (17) 

In case some fraction 77i oss of CR gets lost from the turbulent 
region per cycle, the spectrum evolves towards a power law dis- 
tribution with spectral index 

a= l+^ioss^i < 18 ) 

assuming a continuous injection of low energy CRs is provided. 
For a low loss fraction 77i oss < Xa^, the asymptotic spectrum can 
even become harder than the canonical value a = 2 of diffusive 
shock acceleration in the strongest non-relativistic shocks (and 
in the test-particle limit). 



Given the final duration of the transsonic ICM turbulence 
of 1 Gyr ~ 10 Teddy (partly because the virializing dark matter 
flows continue to stir turbulence for the dynamical time of a clus- 
ter), and the limited pool of low energy CRs, due to the severe 
Coulomb losses in the subrelativistic regime, such a spectrum 
is probably not established in typical galaxy clusters. In a cool 
core, a slightly larger number of eddy turn overs might be realis- 
tic, however, the escape probability of the accelerated CRs into 
the wider ICM might be large, so that the accelerated spectrum 
might be steep. 

This discussion should have made clear that some CR ener- 
getisation can be expected if CR diffusion contributes substan- 
tially to the active CR propagation. In order for this to happen, 
the poorly known CR diffusion coefficient has to be on the right 
order of magnitude to have the diffusion speed being a non- 
negligible fraction of the advection speed. Different to the case 
of CR streaming, we are not aware of any natural explanation for 
this. In case of CR streaming in post-merger clusters, the cluster 
turbulence velocity as well as the CR streaming speed are each 
linked to the sound speed and therefore comparable. 

3.6. Expected CR profiles 

We now generalise the picture to full turbulent advection, com- 
plex magnetic topology, and CR transport. What are the expected 
CR profiles in galaxy clusters that are established by the inter- 
play of advection and propagation of CR? Is the process limited 
to operate only when the assumption of strict CR confinement to 
flux tubes is valid? 

In order to answer these questions, we assume CR diffusion 
to enable the CR to change between magnetic flux tubes and 
thereby find paths to more peripheral regions. Thus the CR prop- 
agation needs not to stop after Lb but can reach even the outskirts 
of galaxy clusters, where the infall of matter onto the cluster be- 
hind the accretion shocks prevents further esc ape 

18 This infall has typical Mach numbers of a few to ten. Thus CR 
would need to stream with highly supersonic velocity against this mat- 
ter stream in order to escape, which does not look very likely. However, 
even if some CR would escape this way, we do not believe that the ef- 
fect would be sufficiently large to be necessarily included in our rough 
modeling. First, as we will see in Sect. 13.71 the periods of low cluster tur- 
bulence seem to be sufficiently long that the spatial CR distribution can 
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Without turbulent advective transport counteracting, we 
would expect a completely homogeneous CR distribution within 
the galaxy cluster to establish itself after some time. Since the 
required perpendicular CR diffusion steps might be slow, the 
macroscopic a st should be smaller than the microscopic one 
o-st of the previous sections. If we model the diffusive transport 
through bottlenecks as a spatial decrement in a s t(/)> where r is 
the (radial) coordinate parallel to the CR gradient, we find 



I r 

L Jo 



dr- 



1 



ajf) 



(19) 



Thus, a few bottlenecks with a st (r) <sc 1, at which the CRs spend 
most of their time with low transport velocity, can reduce the 
macroscopic speed by a significant factor to be well below the 
sound speed. In the following we simply write ar st instead of a st 
for the macroscopic streaming Mach number. 

If the diffusion of the CRs out of magnetic bottlenecks dom- 
inates the transport time, any energy dependence of the effective 
diffusion coefficient imprints itself on the macroscopic stream- 
ing speed. We are not going to model this in the following, how- 
ever, we note, that the energy dependence of the transport can 
lead to spatially varying spectral indices of a CR population. 

The cluster-wide CR profile relaxation time of T re i ax ~ 
^as/iAt ~ 10 r c /vgt ~ 2 Gyr/o-st can be comparable to the cluster 
age (where 7? as * 10 r c is the radius of the accretion shock con- 
fining the CR population within the cluster). Thus, turbulent CR 
advection could be very essential to maintain a centrally peaked 
CRp profile, if we want to explain radio halos with the hadronic 
model, or a centrally peaked low-energy CRe population, for re- 
acceleration models to be operative. CR propagation aims at es- 
tablishing a spatially homogeneous CR distribution, but what is 
the preferred profile of adiabatic advection? 

To work this out, we use the following assumptions. 

1 . The cluster is characterised by a mean pressure profile. For 
the sake of simplicity, we ignore pressure fluctuations due to 
the turbulence. 

2. CR propagation operates on small scales, permitting CR ex- 
change between nearby gas volume elements, but not on 
large scales so that the profile is still determined by the ad- 
vective transport. 

3. Magnetic bottlenecks reduce the macroscopic CR streaming 
speed to a level that it is irrelevant for the global profile. This 
assumption will be alleviated later on. 

Under these conditions turbulent advection will completely 
dominate the CR profile. 

Whenever two volume elements come close at some radius, 
they might exchange CR, tending to establish a constant CR pop- 
ulation in any given radial shell[3 During radial advective trans- 
port from radius r to r' the ICM gas with the enclosed CR is com- 
pressed or expanded by a factor X(r — > r') - (P(r')IP(r)) l l y , 
where P(r) is the pressure profile and y = 5/3. The CR rest-mass 



relax significantly, but not fully. Thus, the bulk of the CR distribution 
never propagates to the cluster boundary. Second, in case high energy 
cosmic rays manage to escape against the cluster accretion flow, many 
more lower energy CR would manage to cross the accretion shock, but 
would be swept back, thereby be compressed and further energetised at 
the shock wave. This implies, that the cluster boundary is more a CR 
energy source, than a CR leak. In order to keep our discussion simple, 
we ignore both, potential CR escape through and CR injection and re- 
acceleration at the cluster accretion shock wave. 

19 Or more precisely, on any surface of constant pressure, which we 
identify from now on with radial spheres for simplicity. 



density g(r) — m J dp f(r, p) (or any other advected quantity) 
therefore wants to establish - under the influence of advection 
alone - a profile according to 



e(r) = go 



P(r) 



0o r](r), 



(20) 



where T](r) = {P{r)/Po) l ^ y is the advective CR target profile. 
This is the only steady state profile which fulfills the requirement 
X(r -> r') = g(r')/g(r) = (P(r')/P(r)) l/y for any r and r' . The 
CR normalisation develops therefore a radial profile 



C{r) = Co 



P(r) 



C (v(r)f«, with 



a = ~ 1.33 - - ■ 1.67. 



(21) 



This profile is strongly peaked at the cluster center, and is very 
much in contrast to the spatially constant profile with j3 cr = 
that CR propagation seeks to establish. 

In reality, CR propagation and advection will both shape the 
CR profile. The ratio of their transport coefficients will deter- 
mine the exact equilibrium shape. We will argue that during 
the transsonic turbulence of a cluster merger, we expect both 
processes to operate on comparable strength and therefore we 
expect a profile which is centrally peaked. A Gyr after a clus- 
ter merger, turbulence decays and CR streaming becomes more 
dominant, leading to a flattened profile implying a negligible CR 
population in the cluster center compared to the turbulent advec- 
tion case. Thus, we expect a large variety of CR profiles being 
present, depending on the turbulent history of the cluster. This 
must have implications on radiative signatures of CRs in clus- 
ters. 

In the following, we want to get some analytical insight into 
the dependence of the CR profile on the turbulence level. All 
we aim for is a rough model, which captures the essential de- 
pendencies, and leaves a more accurate treatment to numerical 
simulations. The CR continuity equation for g in the absence of 
sources and sinks can be written as 



dg 

~dt 



+ V -(vg) = 0, 



(22) 



with v — -L> a d + t>di + u s t the CR transport velocity, which is 
composed by advective (u a{ j), diffusive {v&), and steaming (v st ) 
transport velocities. 

We characterise the passive, advective transport via turbu- 
lence as an additional diffusion process with diffusion coefficient 



(23) 



In a stratified atmosphere the effective and averaged advective 
velocity of this (passive, macroscopic) diffusion is given by 



■Uad = -*tu - V - = -K tu V In - 

g j] 



(24) 



What is a bit unusually here is the appearance of the target den- 
sity profile j](r) - (P(r)/Po) 1 ^ 7 , which ensures that the advective 
diffusion alone tries to establish the steady state CR density pro- 
file we derived above in Eq.|20] Its appearance can be understood 
as follows: 

If turbulent advection is dominating the solution, we should 
find q/qq - T] in the equilibrium configuration, according to 
Eq. [20] Thus the effective gradient in the CR population, which 
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Fig. 4: Left: CR density profiles for y tu = 1,3, 10, 30, and 100 (from bottom to top at small radii) including the same number of 
CRs each. Profiles are normalised to £>(0)|y lii=oo . Also the more narrow gas density profile is shown (thick grey line). Right: CR 
normalisation profiles for the same parameters and the gas density profile. (fi c \ = 0.8, /? as = 10, a = 2.5). 



drives the turbulent transport, has to vanish. This can only be 
achieved by a term proportional to 'V(g/r]). Already by dimen- 
sional reasoning one sees that the prefactor has to be -k (u rj/g. 

The counteracting streaming and diffusive CR propagation 
can also be described by velocities. The streaming velocity has 
magnitude y st and is anti-parallel to the CR space density gradi- 
ent. Active CR diffusion can be described very similar to turbu- 
lent diffusion, just with a constant target space density. Thus we 
have 

V * - (25) 



u st = -v st 



and 



1 



Q 



Vg = -Kdi Vln(e), 



(26) 



with Kdi the macroscopically averaged CR diffusion coefficient, 
which depends on k\\, k x and the magnetic topology in a nontriv- 
ial way. We are not trying to give a concrete value for from 
first principles, however, we assume that diffusive transport is 
subdominant but non-zero, since streaming over distances > Lb 
requires that CR switch their guiding magnetic field lines. If dif- 
fusive CR transport would be significant, it would correspond to 
an anisotropy in the CR phase-space function, which would am- 
plify plasma waves on which the CR would scatter. Shortly, the 
diffusive transport is limited by the streaming regime. 

The CR space density becomes stationary for v — 0, and this 
reads in spherical symmetry with radially outstreaming CRs 



d , (g 

V st = Kta — In - 

or \j] 



d 

+ k& Mg). 
or 



(27) 



This is solved by 



g(r) = g () exp 



dr' 



Qo 



rj(r'f e exp 



Km 



1 



Ktu + fdi 1+KdiMu 



Km + K di 

with 

< 1, and 



(28) 



r* = 



Kt u + K di 



Km 



fftu^tu 



/3 e v st 3/3 e a st 



re- 



in the second line f st , K tu , and at^ were assumed to be spatially 
constant. In the following, we adopt a standard cluster beta- 
profile for the thermal electron density n(r) 



n(f) = n 1 + 



(29) 



Here we assume for simplicity a typical value for /3 C \ ~ 0.8 and a 
constant temperature profile|3 With these simplifying assump- 
tions P(r)/Po = n(r)lno and 



rfif) = 1 + 



(30) 



Thus the CR density profile is a bit flatter than the thermal den- 
sity profile. 

The solution of Eq. [28] is only physical for this profile be- 
tween 



2y 



1+ A 1- 



2r c y 
3 ySd p e r. 



(31) 



since at these radii u st changes sign, which we did not model 
in Eq. [27] The CR profile outside r_ < r < r + is actually non- 
stationary. Here, we just set g(r) - g(r ± ) for r < r_ and r > r+, 
respectively, since the dominating streaming in these regimes 
wants to establish a constant CR volume density, which we take 
as an acceptable approximation to the time-averaged density 
there. These radii are approximately given, in case r_ «c r+, 
by 



3y6 cly 8 e / 

3y6cly8 e 



r 



r. 



amXm Pa 

amXmPcX 
; 

a sl y 



— , and 



(32) 



20 We note that this can be easily generalised to a more realistic declin- 
ing temperature profile in the peripheral cluster regions and emphasize 
that our results do not depend on the exact value of the outer slope of 
the density profile and the external cluster regions beyond R 500 - at least 
for streaming distances that are small compared to the virial radius. 
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A centrally enhanced CR profile exists only if r± are real (see 
Eq.[3Tb. which translates into the necessary condition 

atuXta *tu Tgt 2y 
3a st v st r c r ad 3yS c i 

is fulfilled. It means that in order to develop a centrally en- 
hanced CR profile, turbulence has to be above some critical value 
(comparable to the sound speed), and the CR streaming must 
be inhibited by magnetic topology constrains to a value signif- 
icantly below the sound speed. More specifically, the timescale 
T ad = '"oMu = 3 r tu for turbulent diffusion over the distance r c 
should be significantly shorter than the macroscopic streaming 
timescale r st = r c /y st for the same distance. This implies that 
magnetic bottlenecks are critical in lowering the microscopic 
streaming velocity of CR by some finite factor. In order that the 
macroscopic streaming does not disappear (and we do not have a 
switch off mechanism for halos) some cross-field diffusion must 
be present. Since we expect then also some diffusion along the 
magnetic field lines, we can expect some level of diffusive ac- 
celeration of CRs as discussed in Sect. 13.51 

However, in order to also have a stationary situation for the 
CR spectra, we assume that re-acceleration is negligible, i.e. that 
Kt& ^ x tu and therefore /3 e ~ 1. Since then both dominating 
transport mechanisms, turbulent advection and CR streaming, 
are adiabatic in nature, we can translate the CR density pro- 
file into a profile of the spectral normalisation simply by using 
C{r) = C {Q{r)lQ Y«. 

For our beta-profile cluster this yields 

C(r) = Q,(l + ^j exp(^Ar). with (34) 

3 a + 2 

Pc = -PdPgPa, and y6 cr = — ^— 

within r_ < r < r+ and 

C(r) = C(r ± ) (35) 

for r > r + and r < r_, respectively. 

When we compare the different stages of the same galaxy 
cluster during its evolution during and after a cluster merger, 
we normalise the CR populations to have a constant total CR 
number 

Ncr =4tt f " drr 2 —. (36) 
Jo '« 

This is a conservative assumption, since during a cluster merger 
additional CRs will be injected, so that the CR content of pre-and 
post merger cluster should differ by some margin. The resulting 
CR density profiles are displayed in Fig. [4] 

We will use this toy model in the following to illustrate the 
effect of the CR transport processes on CR signatures in clusters 
of galaxies. It is certainly a very idealised picture in many re- 
spects. For example, the quasi stationary configuration implied 
is probably never reached in real clusters. The finite lifetimes of 
any cluster weather phase probably leads to a constantly chang- 
ing CR profile. Nevertheless, it can be regarded as an extreme 
model, being the opposite extreme to the case of strict CR con- 
finement in finite magnetic structures discussed in the previous 
section. Thus we have two scenarios, which probably bracket the 
physical reality in galaxy clusters: 

1. CR are confined in closed magnetic structures, which due 
to their radial extension, connect regions with some pressure 
ratio and therefore can be characterised by a CR compression 



1 . 

ratio X = (P\/P2) r , as discussed in Sect. 13.41 We refer to this 

as the case of confined CR. 

2. CR are able to stream relative freely (with reduced effective 
streaming velocity) throughout the whole cluster volume, 
and manage to establish the stationary profiles described in 
this Sect. 13.61 We refer to this as the case of mobile CR. 

In the following, we discuss the observational signatures of these 
two extreme cases. But first, we should look at the time-scales 
required to establish those steady-state profiles. 

3.7. Time scales 

The profiles calculated in the previous section are mere illustra- 
tions of how the CR distribution might look like. In reality, the 
CR profiles will always be in an intermediate state between such 
profiles, since the cluster turbulence is never stationary and also 
does not exhibit the same velocities at all radii. In order to see if 
these profiles can be used to guide our intuition, we have to see 
if the real CR distribution is able to change significantly within 
a sufficiently short time. 

Radio halo switch off: We assume that at a certain moment 
CR streaming has become the dominant transport mechanism, 
removing CRs especially from the center. In order to diminish 
or extinguish a radio halo, the central CR density in a cluster 
might fall by a factor of X = 2 or X = 10, respectively. If 
we approximate the central CR distribution as a homogeneous 
sphere with radius r c , this requires an increase of the radius to 
r = X 1 ^ r c = 1.26 r c or 2.15 r c , respectively. Thus, the CR 
have to stream a distance of 0.26 ... 1.15 r c , which takes them 
t = (0.06 . . .0.28) /a st Gyr, respectively, for a sound speed of 
c s = 1000 km s 1 and r c = 250 kpc. Note, that the CR normali- 
sation drops by a larger factor X (Q,+2 ^ 3 » X 15 — 3 . . . 30, respec- 
tively, due to the effect of adiabatic energy losses. 

As we will see in Sect. 14.21 this corresponds to a decrease 
in radio flux in the hadronic model by at least a factor X 1 5 3 = 
3 . . . 1000, depending on the central magnetic field strength and 
the expansion factor. Thus, a switch off of cluster radio halos on 
timescales below a Gyr can be expected even in case that a st < 1 , 
as long it is not too low (a st > 0.06 . . . 0.26). 

In the re-acceleration model, the radio halo switches off with 
the radio halo CRe cooling timescale of about 0.1 Gyr the mo- 
ment the re-acceleration stops to operate, irrespective if there is 
CRe transport from the outside or not. 

Radio halo switch on: Both radio halo models, the hadronic 
and the re-acceleration one, require a seed CR population to 
form a radio halo, either protons or electrons, respectively. 
Again, transport over a core radius is required in order to pump 
CRs into the cluster center. The timescale for this is t — r\ //c tll = 
0.73/(o'tu^ / tu) Gyr, and therefore sufficiently short if the turbu- 
lence is sufficiently transsonic. The initial transport of CR into 
the cluster center during a merger is probably poorly described 
by our turbulent diffusion model. It might be proper to assume 
that the CR just have to be advected by a coherent flow through 
the cluster core, so that t = 2r c /v tu = 0.5/atuGyr. However, 
realistic descriptions should use 3-d numerical hydrodynamical 
simulations and are therefore beyond the scope of this initial pa- 
per on this subject. 

To conclude, the assumption of a switch-on/off mechanism 
for radio halos by CR transport mechanism, which operates on 
timescales below a Gyr, is realistic. The above calculated CR 
profiles can therefore provide us with useful guidance, at least 
for the cluster centers. 
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Fig. 5: Left: Gamma-ray emissivity profiles for the CR distributions in Fig. [4] and X-ray emissivity profile of the ICM in grey. 
Emissivities are normalised to the central emissivity of a cluster with y tu = oo. Right: Total gamma-ray flux due to hadronic CRp 
interactions with the ICM nucleons as a function of y tu = r st /T a( j and for a = 2.25, 2.5, and 2.75 (solid, dashed, and dashed-dotted 
lines, respectively). Normalised to L y for y tu = oo. (B c \ = 0.8, /? as = 10). 



4. Implications for thermal and non-thermal 
emission 

In the following we illustrate the consequences of CR transport 
for the gamma-ray and radio halo emission of clusters, and for 
the heating of cluster cool cores using the two scenarios de- 
scribed above: confined and mobile CRs. 

4.1. Gamma-rays 

The gamma-ray emissivity of a power law CRp spectrum as in 
Eq.|Klis 

/ly = A y C£>gaS, (37) 

where A y is provided in App. [A] and depends on the spectral 
index and the gamma-ray energy window considered, and p gas 
is the gas mass density. The gamma-ray emissivity of a cluster 
with constant metallicity (g gas (r) = fi n e (r), with /i the gas mass 
per electron) is therefore 

dVA y =AnA y n J drr 2 C{r)n e {r). (38) 

If we now assume the case of confined CRp, which expand upon 
vanishing turbulence by a factor X as 2, we find the gamma-ray 
emissivity to be reduced by the factor 

U 

-L =X~ s r * 2~ L5 * 0.35. (39) 

We note that the effect of lowering the target density at the lo- 
cation of the relaxed CRs is exactly canceled by increase in vol- 
ume. Thus, only a moderate reduction of the gamma-ray flux 
has to be expected, which is good news for the gamma-ray de- 
tectability of galaxy clusters. 

In case of mobile CR, the difference in gamma-ray luminos- 
ity between turbulent and quiet clusters is expected to be larger, 
but apparently less than an order of magnitude, as can be seen 
from Fig. [5] 



4.2. Radio halos 
4.2.1. Hadronic model 

Hadronic CRp collisions produce also charged pions, which de- 
cay into relativistic e ± . These cool due to inverse Compton scat- 
tering on CMB photons and synchrotron losses and establish a 
steady state spectrum which is steeper than that of the parent 
protons by one, a e = a + 1. The radio synchrotron emission of 
these electrons has a spectrum with index a v = (a e - l)/2 = 
a /2 w 1 . 15 ... 1 .25. The radio luminosity at frequency v per unit 
frequency interval is 

L v = A v f dVC pggas SB I—) 4 , (40) 

where the abbreviations A v and s Bc are defined in App.lBl These 
quantities are unchanged by CRp transport processes, and can be 
regarded to be constant, in our context. e p h is the energy density 
of the photon field on which the CRe scatter. If this is dominated 
by the CMB, it corresponds to an equivalent field strength of 
Bcmb = 3.27(1 + z) 2 ftG, where z is the redshift. At locations 
where the ICM fields are stronger than this, synchrotron cooling 
dominates, and the radio emissivity depends only weakly on the 
magnetic fields in which the CRp reside. Where B is below B cm b, 
IC scattering is the dominant loss channel, and the radio flux is 
very sensitive to the field strength. 

When the cluster turbulence decays, the CRp occupy mostly 
the weaker magnetised regions of their enclosing field structures, 
in the confined CR scenario. Furthermore, magnetic fields tend 
to become weaker, since the magnetic small scale dynamo main- 
taining the fields is less powerful. We model these effects by 

4 s' a' 2 

s B ^s' B = X£, with X tn = = -B, (41) 

fita < 

the ratio of the turbulence after and during the merger. In case the 
magnetic fields are relatively weak (B <K B C mbX tne cluster radio 
luminosity changes due to the CR profile relaxation according to 

hi ~ x- i ( ff+2 > X tu ^ * X- 3 X m 1125 « 0.06 (42) 

Ly 
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In the approximations, we used a - 2.5, X « 2, and X tu = 2. The 
latter is probably an underestimate of the turbulence ratio, see 
iBattaglia et aT] (12009). If the cluster radio emission is dominated 
by regions with strong magnetic fields (B » B cm b), we find 

y- = X- 4 - ~ X" 1 ' 333 X tu 0125 * 0.36. (43) 

Thus, in the confined CR scenario a radio flux decrement of at 
least a factor three, but more plausibly by a larger margin, can 
be expected. 

In case of a mobile CRp population, the variation of the ra- 
dio flux as a function of turbulence strength is much larger. It 
depends on the magnetic field strength, and especially on how 
much the magnetic fields change in response to a change of the 
turbulence strength. Magnetic field generated by a saturated state 
of a small scale dynamo should have e# oc e tu oc a^ u . However, 
since the magnetic field strength probably requires some time to 
reach the saturation state, we model also weaker dependencies 
by adopting 

= 5 »W /*-}*, (44 ) 
8tt n WO/ 

with Bo = 6 /jG and using 6 b = 0, 1, as well as 2 for no, a moder- 
ate and a strong dependence of the magnetic field strength on the 
turbulence level, respectively. Since we expect the streaming ve- 
locity v st to depend strongly on the turbulence level in an inverse 
fashion, y m = v t ulv st might change much stronger than v tu , and 
therefore we expect 6b < 2. The resulting radio luminosities are 
shown in Fig. [6] 

It is obvious from this discussion that a rapid drop in radio 
luminosity after the turbulent merger phase by one order of mag- 
nitude or more, can easily be achieved; moreover, it is actually 
expected. 

4.2.2. Re-acceleration model 

The presence, luminosity, and spectrum of a radio halo in the re- 
acceleration model depends on two things: the availability of a 
low energy CRe pool and the strength of the wavefield providing 
the CRe acceleration. To calculate the latter is a complex task, 
well beyond the scope of this paper. For magneto sonic waves 
it is calculated in Brunetti & Lazarianl (l2007Ll2010l) . where it is 
shown, that if a significant faction of about 20% of the clus- 
ter pressure is in form of such waves (and their spectrum is 
Kraichnan-like and reaches down to small scales), electrons are 
re-accelerated up to 10 GeV, the radio observable part of the 
CRe spectrum in clusters. It should just be noted, that the genera- 
tion of magnetosonic waves with the Lighthill mechanism scales 
with roughly the second power of the turbulence energy den- 
sity (Light hilll 1951 1 19541) and therefore can be expected to vary 
largely between merging and relaxed clusters. 

Since the radio emitting CRe loose their energy within 0. 1 
Gyr, an arbitrary strong decline of the radio power is possible 
in this model, if no CRp population is present replenishing CRe 
via the hadronic mechanism. In an quiescent phase of the cluster, 
any low energy CRe population in the cluster center is probably 
completely diminished by the severe Coulomb losses (see Fig. 
|2J, and has to be re-established by advection from the outer re- 
gions with less severe Coulomb losses during the next turbulent 
merger. Thus, large variations of the radio lumin osity in the re- 
accele ration model can be expected, as shown by iBrunetti et al.l 
(200§. 
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Fig. 7: CR heating profiles due to CR streaming (upper black 
curves) and Coulomb losses (lower black curves) for the CR 
distributions resulting from different streaming velocities as in- 
dicated at curves, the t urbulent heating pro file (grey dashed line, 
calculated according to EnBlin et al] d2007h ). and the X-ray cool- 
ing profile (grey solid line). The adopted cool core parameters 
are no = 0.05 cirT 3 , kT = 3keV, r c = lOkpc, i» tu = 350 kms -1 , 
B e = 1, B cc = 0.4, R cc = 10 r c , a = 2.5, and a CR to ther- 
mal energy ratio of s CI /e t h - 0.15 at the outer radius of the cool 
core, R cc , and less in the center (down to e cr /eth ~ 0.01 for large 
ar st ). Since we assume the gas motions to be mostly convective, 
we set^ tu = 10 to model an efficient radial transport. This im- 
plies a tu = 0.36 and K tu = 3.6 1 30 cm 2 s -1 . We further have 
rtu = 2.4, 24, and 120 for a st = 0.5, 0.05, and 0.01, respec- 
tively. The total X-ray cooling luminosity of the cool core region 
is about 4.4 10 43 erg s _1 . The heating power of the turbulence is 
about the same, but mostly going into outskirts of the cool core. 
The total CR heating power (streaming plus Coulomb) within 
the R cc = 10 r c is 2.9, 2.4, and 0.8 10 43 erg s -1 for a st =0.5, 0.05, 
and 0.01, respectively. 



4.3. Cool core heating 

The energy lost by CR during streaming heats the ICM, since the 
plasma waves excited by the CR dissipate. The origin of this en- 
ergy is, in absence of CR sources pumping fresh relativistic par- 
ticles into the ICM, the kinetic energy of ICM gas flows. Thus, 
the heating due to CR streaming is not providing an additional 
heat source to the ICM, since the kinetic energy of the gas will 
be dissipated via a turbulent cascade anyway otherwise, but CR 
streaming provides a different channel for this energy dissipa- 
tion, with a different spatial footprint. 

As we show in the following, the heating profile of CR 
streaming is centrally concentrated, and therefore may play a 
role in stabilizing cool core regions of galaxy clusters. In cool 
cores, the coldest gas is at the center and is therefore also the 
densest, leading to the shortest cooling time. In order that no 
cooling catastrophe occurs, as observations indicate, this central 
gas has to be heated preferentially. 

In order to get an idea on the heating profile due to CR 
streaming, we adopt the simple picture for the CR cluster pro- 
file of Sect. l3.6l for a model of an individual cool core. Thus we 
describe the cool core with a beta-profile with small r c ~ 10 kpc. 
However, instead of assuming a fixed number of CRs within the 
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Fig. 6: Left: Radio emissivity profiles for the cluster shown in Fig. [4] assuming the same magnetic field profiles with Bo = 6pG 
and 5b - 0. Emissivities are normalised to the central radio emissivity of a cluster with ~y tu - oo. The X-ray profile is shown 
in grey. Right: Total radio flux due to hadronic CRp interactions with the ICM nucleons as a function of y tu = T st /T a( j and for 
different dependencies of the magnetic field strength on the turbulence level as parametrised by Sb = 0, 1, and 2 (solid, dashed, and 
dashed-dotted lines, respectively). Normalised to L v for y m = 10. (B c i = 0.8, R d>i = 10, Bq = 6 pG, a = 2.5, z = 0) 
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Fig. 8: Central CR heating within r c of the cool core in Fig. [7] 
as a function of the CRp streaming velocity. The CR streaming 
heating power decreases with decreasing streaming velocity. For 
a st ~ 0.1 this trend is weakened by the build up of a centrally en- 
hanced CRp population. For ar st ~ 1 the streaming velocity is so 
large, that the advective work on the CRp population is instan- 
taneously dissipated, and the heat input is therefore independent 
of a st in this regime. The build up of the central CR population 
with decreasing a st can be seen in the Coulomb heating rate. 



cool core region, we assume that the non-cool core region pro- 
vides an environmental CR floor with fixed CR density g m , to be 
dragged into the cool core via the gas motion therein. We still 
model this gas motion again as isotropic turbulence, but keep 
in mind, that in reality it may have more a convective structure, 
since it is probably excited by rising radio bubbles from the cen- 
tral galaxies. These bubbles are one way - if not the dominant 



way - of the central galaxy to provide energetic feedback to the 
cool core preventing catastrophic cooling. CRs escaping from 
such bubbles will add to the heating, as discussed in the refer- 
ences given in footnote [4] We ignore here this contribution, in 
order to see the contribution to cool core heating from the inter- 
play of CR advection and streaming alone. 

The energy deposition of a radially streaming CR with mo- 
mentum pmc is 



dE CI (p) . _ 
dp 

2 p 2 m c 2 

3 Vl + P 



2 2 

p m c 
3 ^JTTp~ 2 



r + 2 



dv st 
dr 



(45) 



If we integrate this over the spectrum given in Eq. [6] (and 
thereby restricting our discussion to the energetically more im- 
portant CRp), while introducing a low momentum cutoff po, to 
be on the safe side, we find the heating power to be 



C(r) m c 2 



a — 2 3 - a 



1 (9?At 

2 dr 



(46) 



with S q (a, b) denoting the incomplete Beta-function, and q = 
(1 + />o( r )) -1 ' Within r_ < r < r+ we can drop the dv^/dr term, 
since there w st = const. Outside of this range, the steady state 
solution requires 



"fad 



d . (g 
Km In - 
d r \ri 



-Km 



drj 3B ci K tu r/r c 



T]d r yr c 1 + r 2 /r| 



since we assumed a constant (time-averaged) CRp density pro- 
file there. 

The low momentum cutof0 depends weakly on the position 
via p (r) = p (n(r)/n ) in . We use po(R cc ) = 1, which roughly 
mimics the Coulom b cooling cutoff the CRp spectrum develops 
(EnBli n et al.ll2007l) . 



21 This cutoff can can, however, be well set to p = in a simplified 
treatment as long as a is sufficiently smaller than 3, which should be the 



18 



EnBlin et al.: Cosmic ray transport in galaxy clusters 



As can be seen in Fig. [7] CRp heating is strongly enhanced 
centrally. This is not only due to the central peak of the CRp 
density, but mostly due to the 1/r factor caused by the streaming 
for r_ < r < r+. But even for completely flat CRp profiles (see 
curve a st = 0.5), where this 1/r-factor does not apply, there is 
a centrally enhanced heating. This results from the Pc/V-work 
done by the gas motion while trying to drag in the CRs, which is 
instantaneously dissipated via streaming. The CRp heating de- 
creases with decreasing streaming speed, since without stream- 
ing no energy can be dissipated through the excitation of plasma 
waves. 

Thus, the CRs collect some fraction of the kinetic energy 
distributed through the cool core, and release it preferentially in 
the cool core center, whereas turbulent cascades just heat each 
mass element with the same rate (in our simplified turbulent cool 
core model): 



etu(r) 



gtuW 
Ttu 



7T L, 



■ £>gas(>) 



(47) 



The parameters of our cool core example in Fig.|7]were cho- 
sen to indicate that CRp heating is a potentially interesting mech- 
anism. We could have picked parameters for which this heating 
would be dominant or insignificant in the center (whereas the 
latter is more easily achieved). It is not clear to us at the moment 
if the self -regulation processes in cool cores drive the system to a 
state where CRp streaming heating is significant or not. This de- 
pends to some degree on the value of a st , which we do not know 
for this environment, as well as on the amount of CRp in the sur- 
rounding of the cool core and other cool core parameters. The 
dependence of the central CRp heating rate on ar st is shown in 
Fig.[8]for the same cool core parameters. For a sufficiently large 
streaming speed (ar st > 0.05) the heating rate due to streaming 
is actually relatively insensitive to this speed and dominates the 
central heating in this particular cool core. 

In order to see under which circumstances CRp streaming 
can be an significant heating mechanism, we have to compare the 
dependencies of the cooling and heating rates with cool core pa- 
rameters. To simplify the discussion, we ignore line cooling (as 
we did in Figs.|7]and[8]l and assume maximally efficient stream- 
ing and thus a spatially flat CRp profile within the cool core re- 
gion. With these simplifications, the profiles of X-ray cooling 
sx(r), of turbulent heating s tu (r), and the heating by CRp stream- 
ing e tu (r), are given by 



s x (r) = -Ao(*T)ii!§(l+i*/i2) 



StuW 



H "e,Q 



(l + r 2 /<f 



(48) 



2 iJ2 



. , , yS c i 1%;^ 3 + r7r; 

Sstin = Pcr,0 



(l+r 2 /r c 2 ) 2 ' 



Here we used the formula for the pressure of a CRp spec- 
trum with a powe r-law momentum distribution as given in 
EnBlin et al.1 (120071) . wrote \x ~ m p for the mean molecular 
weight per electron, and introduced the X-ray cooling constant 



case for the expected a = 2.5. In this case 

C(r) u st m c 2 T 2 (±) C ( r ) v m c 2 
est = — p ~ 2.5 

for r_ < r < r+. The numerical factor in the last expression changes to 
1.2 if a cut off at p = 1 is used, which is typical for CRp under the 
influence of Coulomb losses depopulating the sub-relativistic regime. 



Ao = 5.96 10 24 erg s 1 cm 3 (keV) 1/2 for a metallicity of 0.3 so- 
lar. 

Heating by CRp streaming can therefore become compara- 
ble to the cooling in the center if the CRp pressure is relatively 
high, turbulence is strong, the cool core radius is small, and the 
electron density is low. The ratio of CRp streaming heating and 
X-ray cooling scales as 



g s t(0) -Pcr.O Ax QWTtu kT 
\s x (0)\ Pth,o n e r c 

whereas that with turbulent heating scales as 



e st (0) ^ P^o PccXi 



£tu(0) 



th,o a,. 



(49) 



(50) 



From this we expect that heating by CRp streaming is most 
likely to be of importance in weak cool cores, where the central 
electron density is not that extreme, the cool core radius is small, 
and the temperature is still relatively high. A strong convective 
structure of the turbulence with^u > 1 is also very beneficial for 
heating by CRp streaming, especially to boost it in comparison 
to turbulent heating. Thus, we speculate that CR streaming might 
moderate the initial growth of cool cores, whereas a massive cool 
core is more likely stabilised by the dissipation of stronger tur- 
bulence and other processes not discussed here (dissipation of 
weak shock waves, radiative heating, ...). The cool core param- 
eters of Figs. [7] and [8] were actually chosen to represent such a 
case of a weak cool core with strongly convective turbulence. 

Note, that if CRp diffusion would also be a significant trans- 
port mechanism in cool cores, we expect some level of Fermi I 
CRp acceleration to take place, as discussed in Sect. 13.51 This 
might even allow that a sufficiently strong CRp population builds 
up automatically within the cool core to counterbalance central 
cooling, even if the seed CRp population was not very energetic. 
If such an energetic population could be established this way de- 
pends also on the escape rate of CRp from the cool core, a quan- 
tity we do not know with our limited knowledge on the magnetic 
topology in such environments. One might speculate if the pres- 
ence of radio-mini halos in some cool cores has something to do 
with this possible CR acceleration mechanism. 



5. Discussion 

We argue that streaming is an important CR transport mecha- 
nism in galaxy clusters. It tends, as CR diffusion, to establish a 
spatially flat CR profile, and therefore to drive CRs out of the 
cluster core. This can explain why radio halos are not found in 
every cluster, although simulations indicate that a sufficient num- 
ber of CR protons should have been accumulated in clusters for 
this. 

CR advection, on the other hand, tends to produce centrally 
enhanced CR profiles, as any turbulent mixing process in a strat- 
ified atmosphere tries to establish a constant abundance profile, 
so that the density of the advected quantity follows the density 
of the fluid. 

Thus, CR advection and streaming are counteracting trans- 
port mechanisms. Whenever the former dominates, centrally 
enhanced profiles are established, and whenever streaming is 
more important, a flat profile results. The crucial quantity is the 
advective-to-streaming-velocity ratio, y tu = v tu /v st . 

The advective velocities in galaxy clusters are comparable to 
the sound speed during cluster merger, and less when the cluster 
relaxes after the merger. The streaming velocity is poorly known. 
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On a microscopic scale, it might be of the order of the sound 
speed, or even much larger, if the plasma wave turbulence level is 
low. Microscopically it might be reduced by a large factor from 
its microscopic value, due to magnetic trapping of CR in flux 
tubes and the slow cross field diffusion required to escape. Also 
this topological reduction factor for the streaming speed should 
depend on the level of turbulence that is present in clusters. 

Thus, there are three factors simultaneously increasing y tu = 
Vtu/vst when the cluster turbulence increases: the turbulent ve- 
locities increase, the microscopic streaming speed decreases due 
to larger level of plasma waves, and the macroscopic stream- 
ing speed is further decreased due to a more complex magnetic 
topology. Taken together, the combination of these effects should 
produce a significant variation of y tu between merging and re- 
laxed clusters. 

As a result of this, merging clusters should have a much 
more centrally concentrated CR population than relaxed ones. 
This leads naturally to a bimodality of their gamma-ray and ra- 
dio synchrotron emissivities due to hadronic interactions of CR 
protons. Also in the re-acceleration model of cluster radio ha- 
los these transport processes should be essential, since the re- 
accelerated CR electron populations in the dense cluster centers 
is probably too vulnerable to Coulomb losses, to survive periods 
without significant re-acceleration. Transport of the longer liv- 
ing electrons at the cluster outskirts into the cluster center during 
cluster merger would circumvent this problem. 

Although we did not work his out in detail, it should be noted 
that we expect an energy dependence of the macroscopic CR 
streaming speed, which then should lead to a spatial differenti- 
ation of the spectral index of the CRp population and any sec- 
ondary radio halo emission. Such spectral index variation in the 
radio halo should become especially strong during phases of out- 
streaming CRp, i.e. when a radio halo dies due to the decay of 
the cluster turbulence. 

We have also shown that CR streaming in cluster cool cores 
can help to dissipate the turbulent energy preferentially in the 
cool core center, and thereby potentially help to stabilise cool 
cores against a cooling instability. The CR population in cool 
cores required for this could be either dragged into the cool core 
by turbulent transport or be self-maintained by ongoing Fermi I 
acceleration if CR diffusion is also an important CR propagation 
mode. 

To conclude, we have shown that CR transport mechanisms 
seem to be essential to understand the non-thermal content and 
radiative signatures of clusters. Additionally, they might play a 
role in shaping the clusters' thermal structure, especially in cool 
cores. Further investigations of CR transport in galaxy clusters 
using detailed three-dimensional numerical simulation of the in- 
volved processes are therefore necessary as well as an improved 
understanding of the dependence of the macroscopic CR stream- 
ing velocity on cluster weather conditions. 
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Appendix A: Gamma-ray emissivity 

The gamma-ray emissivity in photon number between the en- 
ergies E\ and Ei of a power law CR p spectrum as in Eq. [6] 
is according to IPfrommer etafl (120081) (which we follow here 



closely) 
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where «n = «h +4«He = Pga.s/»ip is the target density of nucleons 
in a fluid of primordial element composition, 
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the spectrally weighted hadronic cross section, m^oc 2 - 
135 MeV the rest mass of a neutral pion, and the shape parameter 



0.14 or! 



-1.6 



• 0.44. 



(A.5) 



Uy , . ^y 

There is a detailed discussion in IPfrommer & EnBlinl (|2004a) 
how the y-ray spectral index a y relates to the spectral index of 
the parent CRp population a. Here, we assume a y = a. Thus we 
find 



4 m„occr pp 
3 a5 y 



a + 1 a— I 



2 5y 26y 



J /Jje, 



(A.6) 



Appendix B: Hadronic radio-synchrotron emissivity 

The coefficient in Eq. [40] characterizing the radio-synchrotron 
emissivity of electrons injected from the decay of charged pi- 
ons produced in hadr onic CRp-gas collisions is according to 
IPfrommer et al.1 ([2008) (which we follow here closely, but with 
slightly modified normalisation) 

t ^(a e -2)cr TSBc m p \m e ) \GeV J 

with the dimensions [A v ] = erg cm 3 s _1 Hz 1 and the volume 
integral extends over the entire cluster. 

The frequency v dependent characteristic field strength B c , 
and the scaling factor of the emitted energy per time v _1 and 
frequency Ae s nch are given by 



2n ml c 5 v 



327T TOpC 2 a e + 1 



^ 31 (ck)^ (R2) 

(B.3) 



3 e GeW 1 

V^r B c e 3 a e + | F (^) T (^ff) F 



Here V(a) denotes the T-function, a v - (a e - Fj/2 = a/2, and 
B c denotes a (frequency dependent) characteristic magnetic field 
strength which implies a characteristic magnetic energy density 
e Bc =B 2 /(8 7r). 
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